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NOMENCLATURE 


A ,  characteristic  drainage  constant,  m  /  5 

ap,  cross-sectional  area  of  channel  (Plateau  border),  m2 

B ,  characteristic  drainage  constant,  m~  /  5 
Cj,  integration  constant 

C2,  integration  constant 

cv,  parameter  to  account  for  the  mobility  of  channel  walls;  defined  as  the  ratio  of  the  average 
liquid  velocity  in  a  given  channel  with  mobile  walls  to  the  average  liquid  velocity  in  the 
same  channel  if  the  walls  are  immobile,  dimensionless 

Db,  diameter  of  a  bubble,  m 

Ex ,  foam  expansion  ratio,  dimensionless 

G ,  geometrical  constant  equalto  (15  +  7V5)/4  «  7.66 ,  dimensionless 
g ,  gravitational  acceleration,  9.8  m  /  s2 

H ,  height  of  foam  column,  m 

hf,  height  of  moving  boundary  at  the  top  of  the  foam,  m 
hf0,  initial  height  of  the  top  of  the  foam,  m 

hw,  height  of  the  moving  boundary  at  the  top  of  the  liquid  layer  beneath  the  foam,  m 

hw0,  initial  height  of  the  top  of  the  liquid  layer  beneath  the  foam,  m 

L ,  channel  (Plateau  border)  length,  m 

Md ,  mass  of  liquid  drained  from  the  foam,  kg 

np,  number  of  channels  (Plateau  borders)  per  bubble,  # 

P  ,  pressure  of  gas  in  bubble.  Pa 


Plt  pressure  of  liquid  in  a  channel  (Plateau  border).  Pa 

r ,  radius  of  curvature  of  the  interface  between  a  gas  bubble  and  a  liquid  channel,  m 

R ,  radius  of  cylindrical  foam  column,  m 

S ,  cross-sectional  area  of  foam  column,  m2 
t ,  time,  5 

t0 ,  initial  time,  s 

^ ind  ,  induction  time,  5 

tfm ,  time  required  to  fill  a  column  with  foam,  s 

Vb,  volume  of  gas  within  a  bubble,  m3 

Vd ,  volume  of  liquid  drained  from  a  foam,  m3 

v ,  liquid  velocity  within  a  channel  averaged  over  the  cross-section  of  a  channel  (Plateau 
border),  ml  s 

v  ,  liquid  velocity  within  a  channel  averaged  over  the  height  of  the  foam  layer,  mis 
vinj ,  injection  velocity  of  foam,  m  /  s 

vs,  superficial  liquid  velocity,  ml s 

z ,  vertical  coordinate  along  the  height  of  the  foam  column,  m 


GREEK  SYMBOLS 


a  ,  liquid  volume  fraction  of  foam,  dimensionless 

a  ,  liquid  volume  fraction  of  foam  averaged  over  the  height  of  the  foam  layer,  dimensionless 
af,  liquid  volume  fraction  of  foam  injected  into  foam  column,  dimensionless 

aw,  liquid  volume  fraction  at  the  interface  between  the  foam  and  drained  liquid  layers, 
dimensionless 

y,  surface  tension,  D.D225Nlm 


v 


8 ,  model  constant  equal  to  approximately  0.21,  dimensionless 
8a ,  model  constant  equal  to  y/3  -nil  ~  0.161,  dimensionless 
r/ ,  reciprocal  (multiplicative  inverse)  of  a  ,  dimensionless 
6 ,  model  constant  for  comparing  foam  drainage  models,  dimensionless 
p ,  viscosity  of  liquid,  0.000893  kg  /  ms 
permeability  coefficient,  dimensionless 
p ,  mass  density  of  liquid,  997  kg  /  m3 

® ,  generic  function  in  the  standard  form  of  a  linear  ordinary  differential  equation 

y/ ,  generic  integration  factor  used  in  the  analytical  solution  of  a  linear  ordinary  differential 
equation 

Q ,  generic  function  in  the  standard  form  of  linear  ordinary  differential  equation 
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Liquid  Loss  from  Advancing  Aqueous  Foams  with  Very  Low  Water  Content 


1.0  INTRODUCTION 

Filling  a  space  with  foam  is  an  essential  part  of  foam  applications.  One  filling  method  is 
to  apply  foam  directly  onto  a  surface  from  above,  where  fresh  foam  is  added  to  the  top  of  the 
advancing  foam  bed.  This  method  of  filling  is  used,  for  example,  in  fire-suppression 
applications.  To  suppress  fires  in  contained  spaces  such  as  warehouses,  hangar  bays,  and 
obstructed  areas,  aqueous  foams  with  very  low  water  content  are  applied  from  foam 
generators  located  near  the  ceiling,  quickly  filling  the  space.  In  addition,  aqueous  foams  with 
greater  liquid  content  are  sprayed  on  top  of  flammable  liquids  to  suppress  pool  fires  and  build  a 
foam  bed  above  the  pool  to  prevent  re-ignition. 

Another  method  of  filling  a  container  with  foam  is  by  introducing  gas  bubbles  into  a 
liquid  solution  containing  surfactants.  The  gas  bubbles  rise  to  the  surface  of  the  liquid  solution, 
forming  a  foam  bed  above  the  liquid.  In  this  case,  the  fresh  foam  is  added  to  the  bottom  of  the 
column.  Examples  of  applications  employing  this  method  include  foam  fractionation  used  by 
pharmaceutical  and  food  industries  for  protein  separation,  and  froth  flotation  used  by  the 
mining  industry  for  mineral  separation. 

The  loss  of  liquid  from  foam,  also  called  foam  drainage,  has  a  crucial  impact  on  its 
effectiveness  in  applications.  The  dominant  force  responsible  for  drainage  is  gravity,  which  acts 
on  the  foam  from  the  moment  it  is  generated.  However,  most  drainage  models  describe  the 
drainage  of  liquid  from  foam  by  starting  with  a  stationary  foam  column  of  fixed  height.  This 
approach  neglects  the  effect  of  gravity  during  an  essential  stage  of  foam  applications,  the  filling 
process.  Significant  liquid  loss  could  occur  over  the  filling  time,  which  would  ultimately  alter  the 
effectiveness  of  the  foam. 

Modeling  efforts  aim  to  understand  foam  drainage  by  considering  the  flow  of  liquid 
through  the  foam.  The  liquid  flow  in  foam  depends  on  the  geometry  of  the  liquid  contained  in 
the  interstitial  space  between  foam  bubbles.  The  foam  bubbles  take  on  shapes  similar  to 
polyhedra.  In  Figure  1,  an  illustration  depicts  a  bubble  shaped  as  a  dodecahedron  and  the 
surrounding  liquid  structures:  films,  channels,  and  nodes.  The  films  are  thin  sheets  of  liquid 
formed  where  the  surfaces  of  two  bubbles  meet.  The  channels  are  formed  at  the  intersection 
of  three  bubbles:  channels  are  also  called  Plateau  borders  in  honor  of  pioneer  foam  scientist  J. 
A.  F.  Plateau.  The  cross  section  of  a  channel  resembles  a  triangle  with  edges  curved  toward  the 
interior.  The  channels  can  be  considered  as  a  network  of  pipes  that  allow  liquid  to  flow  through 
the  foam.  The  junctions  where  the  channels  intersect  are  called  nodes;  and  four  channels  meet 
at  a  single  node.  The  dimensions  of  films,  channels  ,  and  nodes  depend  on  the  bubble  size  and 
the  liquid  content  of  the  foam.  The  amount  of  liquid  contained  in  the  foam  can  be  quantified  in 
terms  of  the  foam  expansion  ratio  ( Ex)  :  the  ratio  of  the  foam  volume  to  the  volume  of  liquid 

Manuscript  approved  November  8,  2010. 

1 


solution  used  to  form  the  foam;  or  the  liquid  volume  fraction  ( a )  :  the  ratio  of  liquid  volume  to 
total  foam  volume  (Note  that  Ex-\la). 


Figure  1:  Illustration  of  films,  channels,  and  nodes  for  an  idealized  dodecahedral  bubble. 


Most  of  the  modeling  and  experiments  on  foam  drainage  have  been  performed  on 
stationary  foams  with  low  expansion  ratios.  Leonard  and  Lemlich  [1]  modeled  steady-state 
drainage  from  a  foam  fractionation  column,  where  foam  is  continuously  formed  at  the  base  of 
the  column.  The  modeling  approach,  which  is  the  basis  for  most  models,  involved  a  description 
of  the  liquid  flow  within  the  channels  by  considering  forces  due  to  gravity,  variation  in  liquid 
(capillary)  pressure  along  the  height  of  the  column,  and  viscous  resistance  to  liquid  flow.  The 
latter  varied  with  a  parameter  termed  surface  viscosity,  which  was  used  to  quantify  the 
mobility  of  the  channel  walls. 

Desai  and  Kumar  [2]  examined  the  semi-steady-state  liquid  fraction  as  a  function  of 
height  within  a  column  under  similar  conditions,  including  the  loss  of  liquid  from  the  films  with 
time.  The  model  contained  channels  with  a  simplified  cross  section,  an  equilateral  triangle.  A 
parameter  was  introduced  to  account  for  the  mobility  of  the  channel  walls  [3].  The  effects  of 
viscosity,  surface  viscosity,  and  input  gas  flow  rate  on  the  steady-state  liquid-fraction  profile 
were  examined  with  experiments  using  the  surfactants  sodium  lauryl  sulfate  (SDS),  lauryl 
alcohol,  and  Teepol.  The  average,  steady-state  Ex  values  ranged  from  about  10  to  100. 
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Ramani,  Kumar,  and  Gandhi  [4]  proposed  a  model  for  drainage  from  a  static  foam 
column  that  accounted  for  the  time-dependent  drainage  of  films  into  nearly  vertical  and  nearly 
horizontal  Plateau  borders,  where  gravity  was  a  driving  force  only  for  the  liquid  contained  in 
the  nearly  vertical  channels.  The  model  was  compared  with  experiments  using  the  surfactant 
SDS  with  Ex~  8  and  bubble  diameters  less  than  7xlCT4  m  .  Podual,  Kumar,  and  Gandhi  [5] 
later  changed  the  bubble  geometry  from  dodecahedral  to  tetrakaidecahedral  along  with  the 
channel  orientation  of  the  model,  and  the  adjustments  slightly  improved  the  comparison  with 
experiment  for  small  bubble  sizes.  The  experiments  included  foams  with  expansion  ratios  up  to 
100  and  bubble  diameters  up  to  7x10  3  m  . 

The  model  of  Verbist  [6]  contained  a  foam  drainage  equation  which  included  non- 
steady-state  variation  in  the  cross-sectional  area  of  the  foam  channels,  allowing  the  liquid 
fraction  of  the  foam  to  change  with  both  position  and  time.  The  ability  to  calculate  a  time- 
varying  liquid  fraction  including  the  changing  area  of  the  channels  allowed  the  model  to  study 
the  free  drainage  problem:  to  describe  drainage  from  a  fixed-height  column  of  foam  with  a 
liquid  fraction  that  is  initially  uniform.  The  model  viewed  the  foam  channels  as  a  set  of 
independent  vertical  pipes  with  strictly  immobile  walls.  As  in  most  recent  models  of  drainage, 
the  contribution  of  the  films  was  neglected.  Channels  were  considered  in  the  equations  for  the 
conservation  of  liquid  volume  and  momentum  of  the  model,  where  an  effective  viscosity  (about 
150  times  that  of  water)  was  introduced  to  account  for  the  cross-sectional  geometry  and 
random  orientation  of  the  channels. 

The  approach  of  Koehler  [7]  viewed  foam  drainage  as  the  flow  of  liquid  through  a 
porous  medium,  such  as  a  filter.  From  this  perspective,  the  proportionality  between  the  flow 
rate  of  liquid  through  the  foam  and  the  pressure  gradient  that  drives  the  liquid  flow  is 
determined  by  a  permeability  function,  which  depends  on  the  local  liquid  fraction.  The  model 
calculations  were  compared  with  experimental  results  using  the  surfactant  SDS  and  gaseous 
C2F6  to  minimize  the  diffusion  of  gas  between  bubbles.  The  expansion  ratio  studied  for  free 
drainage  from  a  static  foam  column  was  200  and  the  bubble  diameter  was  approximately 
4x10  3  m  .  The  work  suggests,  in  contrast  to  the  model  of  Verbist,  that  the  viscous  effects  could 
be  predominantly  acting  in  the  nodes  rather  than  the  channels  for  foams  with  highly  mobile 
walls. 


Magrabi  et  al  [8,9]  investigated  free  drainage  from  AFFF  and  FFFP  compressed-air  foams 
used  in  fire  suppression.  The  foams  initially  had  an  expansion  ratio  of  20  and  a  bubble 
diameter  of  approximately  2.6x  1(T4  m  .  The  model  contained  a  wall  mobility  parameter  similar 
to  that  of  Desai  and  Kumar  [2,3],  The  work  suggests  that  an  effect  not  included  in  the  model, 
the  increase  in  average  bubble  size  caused  by  gaseous  diffusion  between  bubbles,  could 
significantly  increase  the  drainage  rate  in  these  foams  after  a  few  minutes. 

A  few  recent  works  have  investigated  forced  and  pulsed  drainage.  These  processes 
involve  the  continuous  and  finite  addition  of  liquid  to  the  top  of  a  foam,  respectively.  Drainage 
under  these  conditions  has  been  studied  by  the  non-steady-state  models  of  Verbist  [6]  and 
Koehler  [7],  In  forced  drainage,  the  liquid  volume  fraction  becomes  uniform  behind  a  traveling 
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liquid  front  that  moves  through  the  foam  at  a  constant  velocity.  As  shown  by  the  two  studies 
[6,7],  the  front  velocity  displays  different  power-law  dependence  on  the  liquid  fraction 
depending  on  the  surfactant  solution  used.  However,  forced  and  pulsed  drainage  models  do 
not  apply  to  the  case  of  filling  a  container  with  foam. 

Bhakta  and  Ruckenstein  [10]  have  investigated  filling  a  container  with  foam  from  the 
bottom,  as  in  the  process  of  foam  fractionation.  To  account  for  the  filling  process,  the  model 
included  a  single  moving  boundary,  the  foam  front,  and  considered  the  non-steady-state 
behavior  of  the  liquid  content  within  the  foam.  The  experimental  data  of  Germick  et  al  [11]  that 
included  bubble  diameters  of  about  3x10  3  m  to  4x10  ’  m  were  used  for  comparison.  In  the 
study,  the  effect  of  foam  fill  rate  on  the  drainage  rate  was  not  investigated. 

In  this  work,  we  present  theoretical  analysis  and  experiments  for  drainage  from  a  foam 
as  it  fills  a  container  from  above.  We  introduce  to  the  drainage  problem  an  additional  time 
scale,  the  time  required  to  fill  the  column  with  foam.  This  additional  time  scale  changes  the 
effect  of  drainage  on  the  liquid  fraction  of  the  foam.  Two  moving  boundaries  are  treated  in  the 
problem:  the  air/foam  interface  where  foam  is  added  at  the  top  and  the  foam/liquid  interface 
where  the  foam  layer  meets  the  drained  liquid  that  collects  beneath  it.  The  height  of  the  foam 
bed  is  therefore  controlled  by  both  the  rate  of  foam  injection  and  the  rate  of  liquid  drainage. 
We  also  treat  free  drainage  as  a  special  case  of  the  model,  wherein  there  is  a  single  moving 
boundary  at  the  foam/liquid  interface. 

We  focus  our  investigation  on  one  type  of  aqueous  foam  used  in  practice  for  fire 
suppression  generated  with  a  surfactant  designed  for  achieving  high  expansion  (HiEx).  The 
foams  possess  greater  expansion-ratio  values  and  bubble  sizes  than  those  investigated  in  most 
previous  works.  An  important  aspect  of  this  work  is  that  we  provide  analytical  solutions  for  the 
liquid  drained  from  the  foam  and  the  average  liquid  fraction  of  the  foam  as  functions  of  time. 
The  analytical  solutions  offer  the  ability  to  quickly  predict  drainage  and  to  clearly  understand 
how  the  various  parameters  influence  the  behavior.  Similar  to  previous  works,  the  model 
contains  an  adjustable  parameter  called  a  permeability  coefficient  that  accounts  for  the 
geometry  of  the  channels  as  well  as  the  mobility  of  the  channel  walls.  The  solutions  allow  the 
permeability  coefficient  for  a  given  foam  to  be  calculated  by  fitting  to  experimental  data.  We 
show  good  agreement  with  experimental  results  by  treating  the  permeability  coefficient  as  a 
constant.  We  use  the  model  to  examine  the  effect  of  fill  rate,  expansion  ratio,  column  height, 
and  bubble  size  on  the  drainage  behavior. 


2.0  MODEL 

Filling  a  container  from  above  involves  delivering  fresh  foam  to  the  top  of  a  growing 
foam  bed,  and  the  injection  of  foam  is  ceased  when  the  container  is  full.  Both  during  and  after 
the  filling  period,  liquid  drains  out  of  the  foam  and  forms  a  liquid  layer  beneath  the  foam  bed. 
During  the  process  of  drainage,  bubble  films  empty  a  portion  of  their  liquid  to  the  surrounding 
channels.  Liquid  in  the  channels  moves  downward,  passing  through  nodes  into  lower  channels 
repeatedly  as  it  makes  its  way  out  of  the  bottom  of  the  foam.  As  the  liquid  content  of  the  foam 
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diminishes,  the  films  reduce  in  thickness,  the  nodes  shrink  in  volume,  and  the  cross-sectional 
area  of  the  channels  decreases. 

There  are  several  physical  factors  that  drive  the  flow  of  liquid  in  foam.  For  the  foams  of 
interest,  the  forces  governing  fluid  transport  within  the  foam  are  not  in  an  equilibrium  state 
when  the  foam  is  generated.  For  foam  with  sufficient  liquid  content,  the  gravitational  force  is 
initially  dominant  for  the  liquid  in  the  channels  and  nodes,  driving  liquid  downward  within  the 
foam.  Liquid  flow  can  also  be  driven  in  a  foam  when  there  is  a  spatial  variation  in  capillary 
pressure,  the  pressure  difference  between  the  liquid  and  gas  at  the  curved  surface  of  a  bubble. 
The  capillary  pressure  is  caused  by  the  surface  tension  at  the  interface  of  the  two  phases,  gas 
and  liquid,  within  the  foam  which  acts  to  minimize  the  surface  area  of  the  interface. 


Cross-section  of 
liquid  channel 

Channels  are 
formed  in  the 
interstitial  space 
between  three  gas 
bubbles. 


Expanded  view 
of  curved 
gas/liquid 
interface  at  the 
surface  of  a 
channel 


Liquid 


Horizontal 
component  of 
surface  tension  force 

Surface  tension 
force 


Figure  2:  Surface  tension  force  at  the  gas/liquid  interface  of  a  bubble. 


To  see  how  surface  tension  gives  rise  to  capillary  pressure  in  a  foam,  consider  the 
illustration  of  the  curved  gas/liquid  interface  at  the  surface  of  a  channel  in  Figure  2.  The  surface 
tension  force  is  directed  along  the  surface  and  pulls  at  each  end  of  a  surface  length  element  dl . 
The  horizontal  components  of  the  surface  tension  acting  on  a  surface  length  element  dl 
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cancel,  but  the  vertical  components  do  not.  The  net  force  on  each  curved  interface  is  directed 
toward  the  adjacent  gas  bubble  and  its  magnitude  is  inversely  proportional  to  the  radius  of 
curvature  of  the  surface.  The  surface  tension  force  therefore  increases  the  pressure  of  the  gas 
relative  to  the  surrounding  liquid  by  compressing  the  bubbles  at  the  curved  surfaces  of  the 
channels.  The  Young-Laplace  law 


Pg~P,=L  (1) 

r 

provides  a  relation  between  the  gas  pressure,  P  ,  and  the  liquid  pressure,  P, ,  that  is  valid  at 

S  1 

the  gas/liquid  interface.  Note  that  the  pressure  of  the  liquid  in  the  films  approaches  the 
pressure  of  the  surrounding  gas  bubbles  because  the  films  are  very  flat,  i.e.  is  very  large. 

The  curvature  at  the  surface  of  channels  can  vary  throughout  different  regions  in  a 
foam.  Channels  with  more  liquid  have  a  larger  radius  of  curvature,  and  therefore  have  a  smaller 
capillary  pressure,  i.e.  the  liquid  pressure  gets  closer  to  the  gas  pressure.  In  a  draining  column 
of  foam,  the  liquid  content  at  the  top  of  the  foam  becomes  smaller  than  that  of  the  bottom. 
Hence,  a  gradient  in  the  capillary  pressure  is  established,  which  drives  the  flow  of  liquid.  For  a 
draining  column,  the  flow  due  to  the  gradient  in  capillary  pressure  opposes  the  gravitational 
flow. 


Additionally,  the  flow  of  liquid  through  a  foam  depends  on  the  degree  of  slip  at  the 
channel  boundaries.  Slip  is  a  term  used  to  describe  the  flow  of  liquid  relative  to  its  bounding 
surface;  a  no-slip  condition  in  a  channel  means  that  the  liquid  velocity  at  the  channel  walls  is 
zero.  Slow  or  non-moving  liquid  at  the  walls  causes  a  variation  in  liquid  velocity  (i.e.  stress) 
across  the  channel  cross  section,  giving  rise  to  viscous  drag  that  impedes  flow  within  the 
channel.  The  slip  condition  depends  on  the  mobility  of  the  gas/liquid  interface  at  the  surface  of 
the  bubbles.  If  the  surface  of  the  channel  is  able  to  freely  move  along  with  the  flowing  liquid,  it 
would  cause  plug  flow,  a  uniform  liquid  velocity  profile  across  the  channel  cross  section.  The 
slip  condition  for  foam  channel  flow  lies  within  these  extrema  and  can  result  in  different 
drainage  rates  for  the  foam  as  a  whole.  Evidence  suggests  that  the  mobility  of  the  surface  and 
hence  the  slip  condition  strongly  depends  on  the  surfactant.  In  fact,  studies  have  shown  that 
the  drainage  rate  from  a  foam  can  be  significantly  reduced  by  adding  a  different  surfactant  [12], 

A  suggestion  for  the  physical  origin  of  the  mobility  of  the  channel  walls  is  that  the  liquid 
flowing  nearby  could  shear  surfactant  molecules  from  the  gas/liquid  surface  [7],  If  the  shearing 
creates  a  concentration  gradient  on  the  surface,  then  diffusion  of  surfactant  on  the  surface 
could  oppose  liquid  flow.  However,  if  the  surfactant  in  the  flowing  liquid  can  move  to  the 
surface  much  faster  than  the  rate  of  shear,  then  the  surfactant  can  quickly  replenish  the  surface 
without  establishing  an  appreciable  concentration  gradient.  An  alternative  argument  involves 
the  idea  of  viscosity  at  the  surface,  which  reflects  the  forces  of  attraction  between  surfactant 
molecules.  The  strength  of  the  interaction  could  also  be  so  large  that  the  channel  flow  is  not 
strong  enough  to  shear  the  molecules  apart.  If  one  considers  the  films  as  stationary,  then  the 
channel  walls  would  be  immobile  because  they  are  "anchored"  to  the  films.  Weaker 
interactions  may  allow  the  surfaces  to  shear. 
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Further,  in  a  polydisperse  foam,  bubbles  of  different  sizes  differ  in  their  gas  pressure. 
Small  bubbles  exhibit  greater  pressure  and  Henry's  law  states  that  the  solubility  of  gas  in  the 
surrounding  liquid  will  be  greater.  Neighboring  bubbles  with  lower  pressure  will  have 
surrounding  liquid  with  lower  gas  solubility.  Hence,  a  concentration  gradient  of  dissolved  gas  is 
established  in  the  liquid  phase,  driving  the  dissolved  gas  toward  the  larger  bubble(s).  This 
process  is  called  Ostwald  ripening  and  it  makes  an  appreciable  contribution  to  coarsening,  the 
process  of  increasing  foam  average  bubble  size  over  time.  The  time  scale  associated  with 
coarsening  can  vary  significantly  between  foams  [13];  the  study  of  Magrabi  [8]  on  AFFF  foams 
formed  with  compressed  air  shows  that  bubble  size  begins  to  increase  within  just  a  few 
minutes.  Another  factor  that  contributes  to  coarsening  is  coalescence,  where  film  rupture 
between  two  bubbles  produces  a  single,  larger  bubble. 

In  our  model,  we  set  up  the  problem  such  that  the  foam  is  injected  at  a  constant 
velocity.  Further,  we  consider  one-dimensional  drainage  along  the  z-direction  as  depicted  in 
Figure  3.  The  drainage  problem  is  reduced  to  a  single  dimension  because  the  forces  that 
dominate  liquid  transport  are  in  the  vertical  direction  until  the  foam  approaches  equilibrium. 


Foam  input  (af) 


Figure  3:  Illustration  of  drainage  model  including  moving  boundaries. 
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We  consider  two  moving  boundaries.  The  first  boundary,  z  =  hJt),  is  the  air/foam 
interface  that  moves  towards  the  top  of  the  foam  column  at  the  injection  velocity  v.  . ,  defined 
as  the  foam  injection  rate  per  unit  cross-sectional  area  of  the  container.  The  second  boundary, 
z  =  hw(t),  is  the  foam/liquid  interface,  where  the  solution  drained  from  the  foam  pools  beneath 
it.  The  volume  of  liquid  drained  is  Vd  =  Shw(t) ,  where  S  is  the  cross-sectional  area  of  the  foam 
column.  Once  foam  injection  stops,  the  top  of  the  foam  layer  is  fixed  at  the  height  of  the 
column,  H ,  and  the  foam/liquid  interface  is  the  only  moving  boundary.  This  situation 
corresponds  to  free  drainage  of  a  foam  on  top  of  liquid. 

Our  goal  for  the  model  is  to  predict  the  amount  of  water  drained  from  a  foam  and  the 
change  in  the  foam  expansion  ratio  over  time  in  bench-scale  experiments.  Hence,  for  both 
filling  drainage  and  free  drainage,  the  model  solves  for  the  height  of  the  moving  foam/liquid 
interface,  hw(t ),  and  the  liquid  fraction  of  the  foam  averaged  over  the  length  of  the  foam  layer, 

a(t ) . 


Approximations  are  made  in  the  model  for  the  foam  structure.  We  neglect  the  drainage 
of  liquid  from  the  thin  films,  and  we  therefore  only  consider  the  channels.  For  large  bubbles, 
the  portion  of  liquid  in  the  films  is  negligible  [14].  Further,  the  shape  of  bubbles  in  this  work  has 
been  approximated  as  dodecahedral,  where  the  channels  are  formed  at  the  bubble  edges. 
Although  the  shape  leaves  a  small  amount  of  space  unaccounted  for,  the  work  of  Podual  et  al. 
[5]  indicates  that  the  difference  between  tetrakaidecahedral  and  dodecahedral  geometries  in 
drainage  models  is  negligible  for  large  bubbles.  Also,  we  consider  a  monodisperse  foam,  where 
each  gas  bubble  has  the  same  volume.  In  practice,  foams  contain  a  distribution  of  bubble  sizes. 
Similar  to  the  approach  of  Magrabi  et  al.  [8],  we  assume  that  the  drainage  behavior  of  a 
polydisperse  foam  can  be  captured  by  modeling  the  foam  as  monodisperse,  with  a  bubble  size 
equal  to  the  average  of  the  bubble-size  distribution. 

We  maintain  a  constant  bubble  size  in  the  model  and  therefore  neglect  the  effects  of 
both  Ostwald  ripening  and  coalescence.  For  Ostwald  ripening  to  occur,  there  must  be  a 
pressure  difference  between  bubbles  that  yields  a  concentration  gradient  in  the  dissolved  gas. 
For  larger  bubbles,  the  radii  of  curvature  within  the  foam  will  be  larger,  causing  a  decrease  in 
the  capillary  pressure.  If  the  dispersity  in  the  bubble-size  distribution  is  not  large,  then  the 
pressure  differences  between  bubbles  should  be  small.  Hence,  it  is  likely  that  the  phenomenon 
that  drives  Ostwald  ripening  should  be  less  important  for  large  bubbles.  By  the  reasoning 
above,  the  gas  pressure  is  treated  as  a  constant  in  the  model.  The  liquid  pressure  in  the 
channels  is  assumed  equal  to  the  pressure  of  the  liquid  at  the  bubble/channel  interface,  given 
by  the  Laplace-Young  law. 

As  for  coalescence  in  the  model,  recent  work  [13]  shows  that  the  onset  of  coalescence 
depends  primarily  on  a  critical  liquid  fraction  and  is  independent  of  bubble  size.  The  surfactants 
investigated  in  Ref.  [13]  were  TTAB  and  SDBS;  each  foam  studied  exhibited  a  critical  liquid 
fraction  less  than  0.0008  (Ex  =  1250)  .  The  calculations  of  our  model  indicate  that  this  liquid 
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fraction  is  approached  only  near  the  end  of  the  time  for  which  we  are  interested  in  simulating 
the  drainage  process.  Assuming  that  the  foams  in  our  study  have  a  similar  critical  liquid  fraction 
for  the  onset  of  coalescence,  the  effect  of  coalescence  should  be  negligible  in  our  study  during 
the  drainage  period  of  interest. 

The  mobility  of  the  channel  walls  is  not  explicitly  accounted  for  in  the  model.  We  have 
included  an  adjustable  parameter,  the  permeability  coefficient,  which  accounts  for  the  mobility 
of  the  channel  walls,  random  orientation  of  the  channels,  and  corrections  to  the  geometrical 
approximations  of  the  model.  Note  that  most  drainage  models  include  either  an  adjustable 
parameter  for  wall  mobility  or  assume  that  the  channel  walls  are  immobile. 

From  the  partial  differential  equations  (PDEs)  governing  mass  and  momentum  balance 
of  the  liquid  in  a  foam,  we  seek  ordinary  differential  equations  (ODEs)  with  time  as  the 
independent  variable  to  arrive  at  analytical  solutions  that  describe  the  drainage  process.  By 
approximation,  we  solve  for  the  height  of  the  water  layer  beneath  the  foam,  hw,  and  the 
length-averaged  liquid  fraction  of  the  foam,  a  .  A  detailed  description  of  the  solution  is 
provided  below. 

We  develop  the  model  by  considering  the  partial  differential  equations  for  mass  and 
momentum  balance  in  one  dimension.  Treating  the  liquid  as  incompressible,  the  liquid  volume 
balance  in  the  foam  is  given  by 


da  d 
—  +  — 


dt  dz 


(av)  =  0 


(2) 


where  a  =  a(z,t)  is  the  liquid  fraction  of  the  foam  and  v  =  v(z,t)  is  the  liquid  velocity.  The 
liquid  velocity  of  the  model,  discussed  below,  is  an  effective  vertical  velocity  of  the  liquid 
averaged  over  the  cross-sectional  area  of  a  channel.  The  superficial  liquid  velocity  is  vs  =  av , 
which  is  the  volume  of  liquid  flowing  through  a  unit  of  cross-sectional  area  of  the  foam  column 
per  unit  time.  As  illustrated  in  Figure  5,  the  liquid  in  the  channels  travels  through  the  channel 
cross  section  at  velocity  v,  whereas  the  liquid  exiting  the  foam  (or  a  control  volume)  travels 
through  the  foam-column  cross  section  with  the  velocity  vs .  The  balance  in  (2)  is  written  for  a 
control  volume  of  the  foam;  for  the  column  considered,  the  control  volume  corresponds  to  a 
disc  with  the  same  cross-sectional  area  of  the  foam  column,  but  with  a  small  thickness.  The 
previous  expression  states  that  if  the  superficial  velocity  of  the  liquid  coming  into  the  control 
volume  is  different  from  the  output  superficial  velocity,  there  is  a  corresponding  change  in  the 
liquid  fraction  within  the  control  volume. 
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Driving  forces  for  liquid  flow  in  foam 


Lower  liquid 
pressure 


Capillarity 

Force 


Higher  liquid 
pressure 


Figure  5:  Illustration  of  the  forces  that  drive  the  flow  of  liquid  in  foam. 


Next,  we  consider  some  of  the  geometrical  approximations  of  the  model  that  will  be 
necessary  for  the  momentum  conservation  equation.  The  bubble  shape  in  the  model  is  a 
regular  pentagonal  dodecahedron,  and  this  is  used  to  relate  the  geometrical  parameters  of  the 
bubble  and  its  corresponding  channels  to  the  local  liquid  fraction  of  the  foam.  There  are  np 

channels  per  bubble;  for  a  dodecahedral  bubble,  the  number  of  channels  is  equivalent  to  the 
number  of  edges  that  belong  to  a  single  dodecahedron,  n  — 10 .  Each  channel  has  length  L 

and  area  a  ,  and  the  liquid  volume  contained  in  the  channels  per  bubble  is  therefore  npapL  . 

We  neglect  the  additional  volume  of  liquid  contained  in  the  nodes,  which  should  be  very  small 
for  the  high-expansion  foams  of  interest.  The  volume  of  gas  contained  in  a  dodecahedral 
bubble  is  Vb  =  GL3 ,  where  G  =  (15  +  7v5)/4  «  7.66,  and  the  total  volume  of  gas  and  liquid  per 
bubble  is  the  sum  of  these  quantities. 

We  view  the  foam  column  as  a  space  filled  by  a  collection  of  bubbles  and  their 
respective  channels,  where  the  liquid  fraction  of  an  individual  bubble  unit  is  equivalent  to  the 
liquid  fraction  of  the  foam  as  a  whole.  We  therefore  get  the  expression 
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a 


(3) 


npapL 

nPapL  +  Vb 


■  =  a 


nPaPL 


Vh  l-a 


For  a  high  expansion  ratio,  a  « 1 ,  and  we  therefore  make  the  approximation 


npapL 


=  a . 


(4) 


The  cross-sectional  area  of  the  channel  is  related  to  the  radius  of  curvature  by  ap  =  8ar 2 , 

where  Sa  =  V3  -n  /  2  «  0.161 .  By  substituting  this  relation  and  Vh  =  Gls  into  (4),  we  have  a 
relation  between  the  liquid  fraction  of  the  foam  and  geometrical  parameters  of  the  bubbles: 


a  = 


nP5ar'L 

Gif 


(5) 


where  8  ~  0.21 .  This  approximation  was  made  by  Koehler  [7]  for  tetrakaidecahedral  bubbles. 
Note  that  the  relationship  between  cross-sectional  area  of  the  channel  and  local  liquid  fraction 
is 


G 

a  =  — aL  .  (6) 

np 

To  understand  the  dependence  of  the  channel  area  in  (6),  consider  the  liquid  and  gas  content  of 
a  single  bubble.  If  the  volume  of  the  gas  bubble  remains  constant,  increasing  the  liquid  fraction 
of  the  bubble  by  adding  liquid  to  the  channels  will  cause  the  channels  to  swell,  increasing  the 
channel  cross-sectional  area.  Also,  for  a  bubble  that,  regardless  of  size,  maintains  a  constant 
liquid  fraction,  the  channel  volume  is  proportional  to  apL  but  the  total  volume  of  the  bubble  is 

proportional  to  If .  The  liquid  volume  fraction  is  therefore  proportional  to  ap/  Is .  For  the  liquid 

fraction  to  be  constant  as  the  bubble  size  increases,  a p  must  also  be  proportional  to  L 2 . 

Flence,  the  channel  cross-sectional  area  depends  on  both  the  liquid  fraction  of  the  bubble  and 
the  square  of  the  channel  length. 

Next,  we  seek  an  expression  for  the  momentum  balance  for  a  channel  that  is  oriented 
vertically  within  the  foam.  For  a  given  liquid  and  pressure  gradient  across  a  channel,  the 
velocity  of  flow  through  a  channel  will  depend  on  the  channel  geometry  and  the  slip  condition 
at  the  channel  walls.  Deriving  an  expression  for  the  average  liquid  velocity  in  a  channel  is 
complicated  because  the  cross-section  of  the  channels  has  a  complex  shape  and  we  do  not 
predict  or  measure  the  slip  condition.  Given  the  difficulty  of  the  problem,  we  aim  to  show, 
under  simpler  conditions,  that  the  average  velocity  should  be  proportional  to  the  channel  area. 
In  Appendix  A,  we  derive  an  equation  for  the  average  velocity  through  a  cylindrical  channel  with 
a  circular  cross  section  and  a  no-slip  condition  at  the  boundary.  The  result  is  quoted  here: 
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f  V(aP 
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V 
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(7) 
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The  velocity  in  (7)  is  very  similar  to  the  well-known  Hagen-Poiseuille  equation  for  flow  in  a 
cylindrical  pipe. 


The  same  dependence  of  v  on  ap  is  shown  in  the  expression  derived  by  Desai  and 
Kumar  [2]  for  general  wall  mobility  in  channels  with  the  cross-section  of  an  equilateral  triangle 
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The  parameter  cv  is  a  velocity  coefficient  defined  as  the  ratio  of  the  average  velocity  in  a 
channel  to  the  velocity  of  the  same  channel  if  the  walls  are  strictly  immobile. 

In  both  cases,  the  average  velocity  in  the  channel  can  be  described  by 


%ccL2  f  dP,\ 

~ -  PS-^r- 

P  \  oz ) 


(9) 


where  defined  here  as  a  permeability  coefficient,  accounts  for  bubble  and  channel 
geometry,  slip  condition,  and  the  orientation  of  the  channels.  In  this  work,  £  is  an  adjustable 
parameter  of  the  model. 

We  reiterate  that  the  pressure  in  the  liquid  is  found  from  the  Young-LaPlace  law. 


P,=Pg~L  (10) 

r 

which,  at  a  gas-liquid  interface  with  surface  tension  y ,  relates  the  pressure  of  the  gas  in  the 
bubbles,  P  ,  to  the  pressure  of  the  surrounding  liquid  in  the  channel  with  radius  of  curvature, 
r . 


The  pressure  gradient  dPt  / 8z  is  found  by  substitution  of  (5)  into  (10), 
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(11) 


An  expression  for  the  average  velocity  in  a  channel  as  a  function  of  the  liquid  fraction  then 
becomes 
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(12) 


v  = 


gal} 
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y4§  d 

1  ] 

[«  +  L  0Z 

_4cc  J 

Note  that,  overall,  the  velocity  is  proportional  to  the  area  of  the  channel,  which  is  also 
proportional  to  aL2 .  Meanwhile,  the  pressure  gradient  involves  the  gradient  of  1/ r  where  the 
radius  of  curvature  r  cc  cc  Lyfa  .  The  combined  liquid-fraction  dependence  in  the  pressure 

gradient  term  should  therefore  involve  the  gradient  of  the  radius  of  curvature  (i.e.  yfa  ).  This  is 
shown  by  the  simplification 


a 

a  — 
dz 


(13) 


Substitution  of  the  previous  expression  in  (12)  yields  an  equation  for  the  average  liquid  velocity 
in  a  channel 


v 


where  we  introduce  the  constants 


and 


A  = 


tfpg 
V  f*  J 


'  %  Ly  4s' 

V  V  , 


The  units  for  A  are  m/s  (velocity)  and  for  B  are  m2/s  (area/time). 


(14) 


(15) 


(16) 


For  a  given  foam,  (14)  indicates  that  the  average  liquid  velocity  in  a  channel  is 
determined  by  a  dynamic  competition  between  the  liquid  fraction  and  its  variation  along  the 
height  of  the  column.  Equivalently,  this  is  a  competition  between  the  local  channel  area  and  the 
vertical  gradient  in  the  radius  of  curvature.  For  regions  of  foam  with  uniform  liquid  fraction,  the 
average  velocity  is  proportional  to  the  corresponding  liquid  fraction.  For  cases  wherein  a  foam 
is  filled  very  fast  such  that  the  liquid  fraction  is  approximately  equal  to  the  input,  one  should 
expect  a  constant  superficial  velocity  «  A  a2 .  For  slower  filling,  the  liquid  fraction  will  vary 

along  the  height  of  the  column,  decreasing  at  the  top  and  increasing  near  the  foam/liquid 
interface.  For  this  case,  the  effects  combine  in  different  regions.  At  the  end  of  filling,  the  liquid 
fraction  will  quickly  decrease  with  column  height;  the  dynamic  effects  of  the  liquid  fraction  and 
its  gradient  will  oppose  one  another  throughout  the  height  of  the  foam.  Drainage  ceases  at 
equilibrium  when  these  terms  are  equal  in  magnitude. 
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The  average  liquid  velocity  in  a  channel  also  depends  on  multiple  parameters.  For  any 
case,  it  is  proportional  to  the  permeability  coefficient  and  inversely  proportional  to  the  viscosity 
of  the  liquid.  By  definition,  a  greater  permeability  coefficient  increases  the  velocity  of  the  liquid 
in  the  channels.  The  viscosity  plays  a  role  because  we  are  considering  the  case  where  the  liquid 
at  the  walls  is  moving  slower  than  the  rest  of  the  liquid  in  the  channel.  A  more  viscous  liquid 
will  reduce  the  velocity  difference  between  neighboring  liquid  "particles",  and  therefore  reduce 
the  average  velocity  within  the  channel,  overall.  For  a  given  liquid-fraction  profile,  other 
parameters  affect  the  strength  of  the  terms  in  (14)  on  the  competition  that  determines  the 
velocity.  The  constant  A  is  proportional  to  Z2 ,  reflecting  that  the  term  is  proportional  to  the 
channel  area  ap  oc  aL2 .  Further,  the  effect  of  gravity  is  contained  in  A ,  which  always  acts  to 

increase  the  liquid  velocity  flowing  in  a  vertical  channel.  The  value  of  B  determines  the 
strength  of  the  effect  of  liquid  fraction  variation  with  height  on  average  channel  velocity.  In 
contrast  with  A,  B  increases  linearly  with  channel  length  L,  which  indicates  that  A  will  be 
dominant  for  larger  bubbles  sizes.  An  essential  parameter  in  B  is  the  surface  tension,  y ,  which 
is  responsible  both  for  the  existence  and  strength  of  the  capillarity  effect  for  a  given  variation  in 
the  liquid  fraction  with  height. 

Boundary  conditions  must  be  imposed  on  a  .  We  treat  the  case  for  which  a  is  constant 
at  the  boundaries.  During  filling,  the  top  of  the  foam  contains  the  liquid  fraction  of  the  fresh 
input  foam,  af .  The  choice  for  the  liquid  fraction  at  the  foam/liquid  boundary  is  less  obvious. 

Similar  to  the  approach  of  Desai  and  Kumar  [3],  we  assume  that  the  bubbles  are  close-packed  in 
a  face-centered  cubic  (fee)  arrangement.  Further,  the  shape  of  the  bubbles  is  approximated  as 
spherical,  and  the  liquid  fraction  at  this  boundary  is  assumed  equal  to  the  void  fraction  of  close- 
packed  spheres  in  an  fee  arrangement,  a(z  =  hw)  =  0.26 .  The  boundary  conditions  are  also  used 
during  free  drainage  for  comparison  with  the  experiments  of  this  work.  After  filling,  the  liquid 
volume  fraction  at  the  top  of  the  foam  significantly  decreases  from  af .  Flowever,  we  find  that 

reducing  the  value  of  the  liquid  fraction  at  the  top  boundary  below  af  has  a  negligible  effect 

on  the  predictions  of  the  model.  We  maintain  the  same  boundary  conditions  for  filling  and  free 
drainage  to  avoid  unnecessary  complication  of  the  model  parameters. 


The  drainage  model  in  this  work  is  similar  to  other  drainage  models  in  the  literature  [6- 
10],  and  the  equations  governing  the  conservation  of  mass  and  momentum  within  the  foam 
can  be  expressed  in  the  same  form.  The  general  form  is: 


da  d 

- h  — 


dt  dz 


(«v)  =  0 


v=^{[p^a-[rLe  ]£(>/“), 


(17) 


where  6  =  0.46  except  as  noted  below.  The  values  of  the  permeability  coefficient  E,  from  the 
literature  and  from  this  work  are  provided  in  Table  1.  The  large  values  used  to  model  the 
experiments  of  Magrabi  et  al.  [8,9]  indicate  that  the  foams  in  their  study  drain  much  faster  than 
the  foams  studied  in  the  other  works  listed,  which  is  assumed  to  be  a  result  of  channel  wall 
mobility.  The  other  calculated  and  empirical  values  for  E,  are  for  foams  with  immobile  channel 
walls.  For  randomly  oriented  channels,  the  reported  values  for  £  are  smaller  and  vary  at  most 
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by  about  20%.;  the  largest  contribution  to  the  difference  between  these  values  is  the  bubble 
geometry  of  the  model.  In  this  work,  we  also  calculate  £  for  vertically  oriented  channels  with 
no  slip  at  the  channel  walls,  which  is  the  largest  physically  reasonable  value  of  £  for  this  slip 
condition. 


Table  1:  Permeability  coefficients  of  foam  drainage  models. 


Work 

Calculation  method 

£ 

Surfactant  (type) 

Bhakta  and 

Ruckenstein 
(Ref.  [10])* 

Calculated  for  randomly  oriented 
channels  with  triangular  cross 
section  and  immobile  walls0 

0.0055 

Model  results  compared  with 
experimental  data  of  Ref.  [11] 
using  Bovine  Serum  Albumin, 

BSA  (protein) 

Koehler  et  al. 

(Ref.  [7])  t 

Calculated  for  randomly  oriented 
channels  with  scalloped-triangular 
cross  section  and  immobile  channel 
walls1 

0.0063 

N/A 

Magrabi  et  al. 

(Ref.  [8])* 

Empirically  determined  from 
drainage  theory  for  randomly 
oriented  channels  with  triangular 
cross  section  and  mobile  walls0 

0.0245* 

Aqueous  film-forming  foam, 
AFFF:  contains  perfluoroalkyl 
sulfonate  salts  (fluorocarbon 
based) 

Magrabi  et  al. 

(Ref.  [9])* 

Empirically  determined  from 
drainage  theory  for  randomly 
oriented  channels  with  triangular 
cross  section  and  mobile  walls0 

0.0105* 

Film-forming  fluoroprotein 
foam,  FFFP:  contains  hydrolized 
proteins  and  fluorosurfactants 
(fluoroprotein) 

Verbist  et  al. 

(Ref.  [6])0 

Calculated  for  randomly  oriented 
channels  with  scalloped-triangular 
cross  section  and  immobile  walls 

0.0051 

N/A 

This  work 

Empirically  determined  from 
drainage  theory  for  randomly 
oriented  channels  with  scalloped- 
triangular  cross  section  and 
immobile  walls  (Ref.  [7])° 

0.0051 

Buckeye  HXF  2.2%  High- 
Expansion  Foam  Concentrate: 
contains  sodium  laureth  sulfate 
(hydrocarbon  based) 

This  work 

Calculated  for  vertically  oriented 
cylindrical  and  equilateral  triangular 
channels  with  immobile  walls  (see 
Appendix  A)° 

Cylindrical: 

0.03 

Triangular: 

0.022 

N/A 

*Mathematical  form  of  Eq.  (17)  obtained  when  a  «1 


+  0  =  0.41 

t  %  for  initial  liquid  volume  fraction  0.05 

OEquation  (6)  of  this  work  used  to  relate  a  to  a 

°Assumed  dodecahedral  bubble  geometry,  L=0.816R 
TAssumed  tetrakaidecahedral  bubble  geometry,  L=0.72R 
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The  next  step  in  the  development  of  the  model  is  to  obtain  ordinary  differential 
equations  from  the  partial  differential  equations  that  describe  the  balance  of  liquid  volume  and 
momentum  in  the  foam.  To  eliminate  the  spatial  dependence  of  the  equation  for  liquid  volume 
balance,  we  integrate  (2)  with  respect  to  the  vertical  coordinate  z  from  the  height  of  the  liquid 
beneath  the  foam  ,hw,  to  the  height  at  the  top  of  the  foam,  hf : 

hf  Qa 

j  —  dz  +  a(hf,t)v(hf,t)~a(hw,t)v(hw,t)  =  0.  (18) 

K  ™ 


Using  Leibniz's  rule  and  the  boundary  conditions  v(hf,t )  =  v(hw,t )  =  0 ,  this  expression  yields 

h 

—  \adz  -  a(hf,t) — —  +  a(h,t)—^L  =  0.  (19) 

dt  l  v  f  ’  dt  v  H  7  dt 

Let  us  define  a  by  the  expression 


a  = 


(20) 


where  a  represents  the  length-averaged  liquid  volume  fraction.  The  quantity  a  is  useful  for 
providing  an  overall  description  of  the  liquid  quantity  in  the  foam.  Substitution  of  (20)  into  (19) 
yields 


d_ 

dt 


a(hf-hw) 


=  a 


dhf 

dt 


-a 


(M) 


(21) 


where  a(hf,t )  is  the  liquid  fraction  at  the  top  of  the  foam  and  a(hw,t )  is  the  liquid  fraction  at 

the  bottom.  For  the  present  model,  the  liquid  fractions  at  the  boundaries  will  be  considered  as 
constants:  af  =  a(hf,t)  and  aw  =  a(hw,t)  ■  We  consider  that  liquid  drains  from  the  foam  into  a 

liquid  layer  beneath  the  foam,  and  therefore  we  set  aw  =  1  for  the  liquid  volume  balance.  The 
left-hand  side  of  Eq.  (21)  is  the  time  derivative  of  the  liquid  content  within  the  foam,  which  is 
controlled  by  the  quantities  on  the  right-hand  side.  If  we  consider  free  drainage  only, 
dhf  /  dt  =  0  and  the  rate  of  water  loss  can  be  determined  by  the  rise  of  the  water  layer  beneath 

the  foam.  However,  by  considering  the  case  of  filling,  dhf  /  dt  *0,  the  liquid  content  in  the 

foam  is  now  determined  by  the  competition  between  filling  and  drainage.  This  competition  can 
be  controlled  by  the  time  scale  introduced  in  the  model,  the  fill  time. 

Integrating  the  momentum  equation  (14),  we  obtain 
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(22) 


hf  hf 

jv  dz  =  A  ja  dz  -  B  ~ 

hw  hw 

We  define  the  liquid  velocity  averaged  over  the  height  of  the  foam  layer,  V ,  with  the  following 
expression 


hf 

Jv  dz 
—  _  K 

(*/-*.) ' 

Inserting  definitions  (20)  and  (23)  into  (22)  gives  the  expression 

v  (hf  - hw)  =  Aa(hf  - hw ) -  -■&)• 

Note  that  the  boundary  condition  for  aw  in  the  momentum  balance  is  not  the  same  as  that  for 
the  liquid  volume  balance.  The  next  step  is  to  use  these  length-averaged  quantities  to 
determine  an  expression  for  the  height  of  the  water  layer  as  a  function  of  time,  hw(t ) .  The  time 
derivative  of  hw  is  the  volumetric  flow  rate  of  liquid  from  the  foam  layer  into  the  liquid  layer 
per  unit  area  of  the  foam  column.  Equivalently,  this  is  the  superficial  liquid  velocity  exiting  the 
foam. 


(23) 

(24) 


dK 

dt 


■  =  v„ 


(25) 


The  liquid  velocity  and  the  superficial  liquid  velocity,  both  space-  and  time-dependent,  are 
related  by  vs  =av.  If  the  foam  bed  were  separated  into  control  volumes  along  the  height  of  the 
column,  the  superficial  velocity  of  the  liquid  exiting  the  foam  due  to  drainage  would  be 
determined  by  the  liquid  balance  of  the  control  volume  at  the  bottom  of  the  foam  near  the 
foam/liquid  interface.  The  superficial  velocity  of  the  liquid  entering  this  control  volume  will 
depend  on  the  balance  over  the  control  volume  above,  and  the  same  is  true  for  each  control 
volume  in  the  foam.  Our  approach  does  not  resolve  the  foam  column  into  control  volumes; 
instead,  we  seek  a  description  of  foam  drainage  based  on  the  average  properties  of  the  foam. 
We  therefore  make  the  approximation  on  the  superficial  velocity 


dh 


w 


dt 


=  av 


(26) 


and  we  assume  that  the  error  introduced  in  averaging  is  small.  Solving  for  V  in  (24)  and 
substituting  in  the  previous  expression  yields  an  equation  for  the  rate  at  which  the  height  of  the 
water  layer  increases 


dh 


w 


dt 


Aa 2  - 


Ba 


\K-K_ 


(27) 
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The  expression  above  is  used  for  both  filling  and  free  drainage  conditions. 

There  are  two  ordinary  differential  equations  in  the  model  calculations,  one  for  liquid 
volume  balance,  Eq.  (21),  and  one  for  momentum  balance,  Eq.  (27),  where  hw  and  a  are 
unknown.  In  the  following  sections,  two  cases  are  considered  for  the  model:  the  general  case  of 
filling  and  the  special  case  of  free  drainage.  For  each,  ordinary  differential  equations  are 
presented  that  can  be  solved  numerically.  Additionally,  approximations  are  made  to  obtain 
analytical  solutions,  and  the  conditions  for  which  they  are  valid  are  discussed. 

FILLING 


For  filling,  we  integrate  Eq.  (21)  with  respect  to  time  to  obtain  the  expression 

a(hf  - hj  =  afhf  - hw  +  C,  (28) 

where  the  last  term  is  an  integration  constant  representing  the  initial  conditions, 

Ci  —  a(t0)(hf0  —  hw 0)  —  (Xj-hj 0  +  hw0 .  (29) 

The  initial  height  of  the  water  layer  is  set  such  that  hw0  =  0  .  To  prevent  a  singularity  in  Eq.  (27), 

the  problem  is  set  up  at  t  -  0  such  that  there  is  a  small  amount  of  foam  with  height 

hf0  =  0.01  m  having  the  same  liquid  fraction  a(t0 )  =  af  as  the  foam  that  will  fill  the  column.  By 

setting  this  condition,  we  sacrifice  the  description  of  drainage  behavior  within  the  first  few 
seconds  of  foam  application,  but  this  is  not  of  interest  for  this  study.  Under  these  initial 
conditions,  Cx  =  0  .  The  liquid  balance  in  (28)  states  that  the  liquid  content  in  the  foam  is  equal 
to  the  difference  between  the  amount  of  liquid  added  to  the  foam  and  the  amount  of  liquid 
that  has  drained  out. 

Solving  for  a  in  (28)  yields 


_  aih,  ~K 
a  = 


hf~K 


(30) 


Substitution  of  this  expression  in  Eq.  (27)  provides  the  filling  equation 

dh 


dt 


=  A 


afhf-K 

hf~K 


-B 


arhf-K 

_hf~K\ 


(2  -ja, 


(31) 


The  previous  expression  is  used  for  the  numerical  solutions  during  filling  in  this  work. 
We  also  provide  an  analytical  solution  below.  The  numerical  solutions  were  determined  by  the 
function  NDSolve  (MATHEMATICA,  Wolfram  Research)  with  default  settings.  NDSolve  utilizes  a 
variant  of  the  Livermore  Solver  for  Ordinary  Differential  Equations,  LSODA,  to  numerically  solve 
ordinary  differential  equations.  Both  stiff  and  non-stiff  systems  can  be  detected  and  solved 
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using  LSODA,  where  backwards  differentiation  formula  (BDF)  methods  and  Adams  methods  are 
used,  respectively. 

Approximations  are  necessary  in  order  to  obtain  an  analytical  solution  for  filling.  First, 
we  perform  a  Taylor  expansion  of  a 2  about  the  value  a  f 

a2  »a2 \+2af^dc-af~^  =  2afa-a2f  (32) 


The  value  of  a 2  must  be  greater  than  zero,  so  the  approximation  clearly  fails  if  the  condition 
a  /  af>  0.5  is  not  satisfied.  Note  that  this  expansion  is  equivalent  to  neglecting  the  term 

proportional  to  (hw  /  hf  )2  in  Eq.  (31).  Physically,  the  approximation  is  valid  when  a  foam  drains 
much  less  than  50%  of  its  liquid  content  as  it  fills.  We  substitute  Eq.  (30)  into  Eq.  (32)  to  obtain 


a 2  «  2  af 


afhf  ~  K 

hf~K 


(33) 


and  substitute  Eq.  (33)  into  Eq.  (27)  to  obtain 


dK 

dt 


=  A 


2a, 


QCfhf-h  w 

hf~K 


-a, 


-B 


afhf~K 


a,., 


(34) 


To  obtain  a  linear  equation,  we  make  the  approximation  hf»hw.  The  maximum  value  of  hw  is 
afhf,  so  the  approximation  is  equivalent  to  af«  1.  This  approximation  is  valid  for  high- 
expansion  foams,  where  af  is  typically  0.01  or  less.  Using  this  approximation,  Eq.  (34)  can  be 
expressed  as 


dh 


w 


dt 


A 


or, 


2arK 

hf 


(35) 


Putting  Eq.  (35)  in  standard  form  yields 


dK+ 

dt  h. 


Aa2 


(36) 


a  linear  ODE  which  can  be  solved  analytically.  For  mathematical  simplicity,  we  make  one  more 


approximation  in  the  linear  coefficient;  we  neglect  the  term  ( B 


a,.,  -Jcr, 


/hi) .  Since 


hf  =  hf0  +  vinit  is  linear  with  respect  to  time,  the  term  2 Aafhf  1  is  inversely  proportional  to  t, 


‘/o 


/"/ 


whereas  B 


a ...  -  Aa, 


h2  is  inversely  proportional  to  t2 .  Flence,  the  term 
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B  yfa^-ylaj  hf 2  quickly  becomes  much  smaller  than  2 Aafhf  1 .  The  ordinary  differential 
equation  from  which  we  obtain  an  analytical  solution  is  therefore 


Eq.  (37)  can  be  re-written  as 


Eq.  (38)  indicates  that  the  drainage  rate  will  be  equal  to  the  difference  between  the  drainage 
rate  of  a  foam  that  has  a  =  af  (see  Eqs.  (26)  and  (24))  and  a  correction  factor  that  depends  on 

the  "volume  per  cross  section"  ratio  of  the  liquid  lost  to  the  liquid  layer  hw  to  the  total  liquid 
added  to  the  foam  afhf .  The  correction  factor  accounts  for  the  fact  that  a  decreases  from  af 

as  liquid  is  lost  from  the  foam  to  the  liquid  layer  beneath.  Intuitively,  the  ratio  of  the  volume  of 
liquid  in  the  liquid  layer  to  the  volume  of  liquid  in  the  foam  layer  will  depend  on  how  the  rate  of 
drainage  compares  with  the  rate  by  which  the  column  is  filling  with  foam.  As  shown  below,  the 
competition  between  these  rates  is  captured  in  the  analytical  solution. 

Equation  (37)  is  in  the  form  dhw  /  dt  +  ®(t)hw  =  Cl(t).  Using  the  integration  factor 
y/  =  expQ the  solution  is  given  up  to  a  constant  by  hw  =  y/~l  J Q(t)ys  dt .  The 
integration  factor,  using  dhf  =  vinjdt  ,  is 


2  Acif 

=  h  Vlni 

nf 


(39) 


The  solution  is  given  by 
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v  f  ?  -fidf  f  r  I -  I - n  v 

hw=hf  inj  Aaf - —  Jaw-Jaf  hf  inj  dt  +  C ' 


-2Aar  2  2 Aa^ 

v  .  Act  r  r  v.  . 

;n;  2  I  JL  inJ 


2  Aaf  rt  I  /  2.4a/- 

.  — L  Baf  Jaw-Jaf  r  — L-\ 

J  hf  v«  dhf - CL1 - hf  'w  dhf  +  c2hf 


Aa}  vmjhf  % 
vinj  2  Aaf  +  vinJ 


-hf  v* 


BaA  Ja„,  -  Aa 


f  LV“w  V“/J  Vinjnf 


_  V*A 


+  C,A, 


2Aaf+vinj 


+  C2hf 


i  ,  5[V^-V«7]  : 

-  arhf - - - -  +  C2nf 

v  ’  ’  ?  a  1 

2  +  ^l_ 


We  use  the  initial  condition  hw{  0)  =  hw0  =  0  and  hf(  0)  =  hf0  at  time  t  =  0,  which  yields 


0=  -  afhrn - - - -  +  C,A, 

„  v„,  /  /0  2/4  2  f 

2  +  —^— 

Aaf 


which  can  be  used  to  solve  for  the  integration  constant 


afhf0  h/o 


2h — — 
Aa r 


Substitution  of  (42)  into  the  final  expression  in  (40)  yields  the  analytical  solution 
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h  = 


v. . 

2  + 

inj 

Aaf 

1 

2  + 

Vinj 

B 


afhf  -- 


law-Jaf 


2  A 


B 


K~^ai 


2  A 


v. . 

2  +  -  m] 


Aa, 


~ 

2  Actf 

afhf0 

h  Vinj 

70 

.  . 

-2  Aaf 


Aa, 


«/V+ 


v. . 

2  +  -  "y 


Aa, 


B 


afhfo 


K~Aar 


2  A 


(43) 


1- 


lAcij 

hf0 

Vinj 

1 

o 

+ 

The  value  of  a  during  filling  is  found  by  substitution  of  Eq.  (43)  into  Eq.  (30). 


Equation  (43)  can  be  simplified  by  neglecting  hf0 .  The  abbreviated  analytical  solution 
for  filling  is  given  by 


1 


v.  . 

2  +  ^- 
Aaf 


«/V“ 


(44) 


For  cases  where  the  analytical  solution  is  valid,  the  abbreviated  solution  given  in  Eq.  (44)  shows 
very  little  difference  with  the  solution  of  (43)  after  roughly  10  seconds.  The  first  term  in  the 
previous  expression  quickly  becomes  dominant  during  the  filling  process.  This  indicates  that  hw 
has  a  predominantly  linear  time  dependence,  and  therefore  the  drainage  rate  from  the  foam  is 
approximately  constant  during  filling.  The  first  term  in  Eq.  (44)  contains  vinJt ,  which  is  the 

height  of  foam  added  from  the  start  of  filling  to  the  time  t .  The  liquid  portion  of  the  added 
foam  is  afvinjt .  The  coefficient 


1 


v.  . 

2-t — — 

Aaf 


(45) 


is  therefore  a  percentage  of  the  added  liquid  that  is  drained  (and  added  to  hw )  at  time  t .  For  a 
given  input  liquid  fraction  a,,  the  key  parameters  that  determine  this  percentage  are  v  .  and 

J  '■•V 

A  ,  which  represent  the  time  scales  for  filling  and  drainage,  respectively.  A  key  point  of  this 
study  is  that  an  appreciable  amount  of  the  added  liquid  can  drain  when  the  filling  and  drainage 
time  scales  are  comparable.  However,  as  one  might  expect,  very  little  of  the  added  liquid  is 
drained  during  filling  if  one  fills  very  fast  such  that  v.  .  »  Aa f  .  The  analytical  solution  is  not 

valid  for  very  slow  foam  injection  such  that  Aaf  »  vinj  because  the  condition  for  the  validity  of 
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the  Taylor  expansion  of  a 2  is  not  satisfied,  i.e.  a  will  not  be  much  less  than  50%  of  af . 
However,  the  case  can  still  be  described  by  numerical  solution  of  Eq.  (31). 


FREE  DRAINAGE 

For  free  drainage,  we  set  dhf  /  dt  =  0  in  Eq.  (21)  and  integrate  with  respect  to  time  to  obtain 

a(H-hw)  =  -hw  +  Cl,  (46) 

where  we  solve  for  the  integration  constant 

C,  (47) 

Note  that  free  drainage  occurs  from  a  column  of  fixed  height,  H  .  Eqns.  (46)  and  (47)  simply 
state  that  the  liquid  content  in  the  foam  is  equal  to  the  difference  between  the  initial  liquid 
content  and  the  amount  of  liquid  that  has  drained  out,  causing  the  liquid  layer  to  rise  by  the 
height  hw-hw0. 

We  consider  the  general  case  for  initial  conditions  on  free  drainage.  Solving  for  a  in  Eq.  (46) 
yields 


-  +  ~k*\ 

H-h  , 


a  = 


(48) 


Substituting  the  previous  expression  into  Eq.  (27)  yields  the  free-drainage  equation 


dK 

dt 


=  A 


-B 


hw0-hw  +  a(t0)[H-hw0 

h-K 

Ko-K+a(t0)[H  -Kn\ 


[H-K] 


(49) 


a,„  -Ja, 


Eq.  (49)  can  be  solved  numerically.  However,  we  calculate  an  analytical  solution  that  is  nearly 
identical  to  the  numerical  solution  for  all  cases  examined  in  this  study. 

For  free  drainage,  we  change  the  dependent  variable  of  the  ordinary  differential 
equation  to  a  .  Making  this  change  allows  the  calculation  of  an  analytical  solution  that  is  both 
simple  to  obtain  and  valid  for  a  greater  range  of  parameters  than  the  solution  obtained  using 
hw  as  the  dependent  variable.  Differentiating  Eq.  (46)  with  respect  to  time  yields 


da 

dt 


dh 


dt 


(50) 
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We  substitute  the  previous  expression  in  Eq.  (27)  to  obtain 


da 

dt 


Aa 2  Ba 

H-hw  +  [H-hwf 


(51) 


As  in  the  filling  case,  we  make  the  approximation  that  the  height  of  the  water  layer  is  much  less 
than  the  height  of  the  foam  column,  h  «  H ,  and  divide  by  a 2  in  Eq.  (51)  to  yield 


1  da  _  A 
H 


a 2  dt 


B 


Using  the  relation 


d_ 

dt 


a 


H2 


1  da 
a  2  dt 


1 

a 


(52) 


(53) 


and  defining  77  =  Ha,  we  can  express  Eq.  (52)  as 


dr] 

dt 


B 


a,.,  -Aat 


H 


n  = 


A 

H 


(54) 


Eq.  (54)  is  in  the  form  dr]  I  dt  +  <&(t)r/  =  Cl(t)  .  Using  the  integration  factor  y/  =  expQ ® (t)dt^j, 
the  solution  is  given  up  to  a  constant  by  r]  =  y/  1  dt .  The  integration  factor  is 


yr  =  exp 


B  A la,. ,  -Jat 


H2 


-dt 


V 


=  exp 
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a,.,  a, 


H 2 
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(55) 


and  the  solution  is  given  by 


i]  =  exp 


-B 


law  ~\af 


H 2 
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J  77exp 


AH 
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law-Jaf 


+  C2  exp 
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-B 
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law  -  yjaf 


H2 


dt  +  C-, 
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l°cw-Jaf 
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V  J 

Using  the  initial  condition  r/(t0 )  =  l/«"(70)  at  time  t0 , 


(56) 
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tiM  =  — 
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AH 


aw  ~Jaf 


«Oo)  B 

yields  the  integration  constant 

r  1  ah 

c2 


+  C2  exp 


-5 


locH,-Jaf 


H2 


(57) 


law-Jaf 


<*M  B 

The  analytical  solution  for  free  drainage  is  therefore 


exp 
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law-^af 


H2 


(58) 
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AH 
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exp 
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H1 


(59) 


and  the  corresponding  value  of  hw  can  be  found  by  solving  for  a  in  Eq.  (59)  and  substituting  in 
Eq.  (48). 


Let  us  examine  a  limiting  case  of  Eq.  (59).  As  t  approaches  infinity,  we  have 


1  _  B 
=  a(oo)  =  — 


a -Ja, 


>7  (go) 


AH 


(60) 


The  value  of  a( oo)  is  the  equilibrium  liquid  fraction  of  the  foam  as  predicted  by  the  model.  This 
expression  for  a( oo)  could  also  be  arrived  at  from  Eq.  (27)  by  using  the  condition  dhw  /  dt  =  0  at 
equilibrium  and  the  approximation  hf  -  H  »hw. 


We  note  that  the  average  expansion  ratio  of  the  foam  is  given  by  Ex-\t  a  .  We  can 
express  Eq. (59)  as 


f 

Ex(t)  -  Ex  (oo)  =  ^Ex(t0)-  Ex  (oo)  j  exp 

v 
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<af 
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(61) 


which  shows  that  the  foam  progresses  from  its  initial  expansion  ratio  to  that  of  equilibrium  on 
an  exponential  time  scale.  Note  that  the  exponent  controlling  the  rate  at  which  the  initial 
expansion  ratio  progresses  to  equilibrium  strongly  depends  on  the  height  of  the  column;  the 
process  is  slower  for  taller  columns. 
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3.0  EXPERIMENTS 


To  generate  foam,  we  created  a  liquid  foam  solution  by  mixing  2.2  parts  concentrated 
liquid  foam  solution  (Buckeye  HXF  2.2%  High-Expansion  Foam  Concentrate)  with  97.8  parts 
distilled  water.  The  surfactant  solution  is,  by  weight,  26%  water,  29%  hexylene  glycol,  38% 
sodium  laureth  sulfate,  and  8%  proprietary  mix  of  surfactants/stabilizers  (refer  to  MSDS).  We 
poured  the  solution  into  a  liquid  chamber  and  pressurized  the  chamber  with  nitrogen  gas  ,  N2. 
The  liquid  solution  was  fed  to  the  foam-generation  device  illustrated  in  Figure  6.  The  liquid 
solution  entered  the  side  of  the  foam  generator  and  exited  downward  through  a  cone- 
patterned  nozzle  (Model  302-C,  SureShotsSprayer.com)  onto  a  copper  screen  with  a  diameter 
of  6.4  cm  and  30x30  mesh  cells  per  inch  (40.8%  open  area)  held  in  place  by  a  rubber  gasket.  The 
N2  pressure  was  increased  in  the  liquid  chamber  such  that  the  nozzle  sprayed  the  liquid  solution 
at  a  fixed  volumetric  rate  of  about  60  mL/minute,  measured  with  a  calibrated  rotameter  (Dwyer 
VFA-33  Visi-Float  Flowmeter).  An  in-house  source  provided  air  flow  into  the  top  of  the  foam 
generator.  A  perforated  plate  and  a  bed  of  tightly  packed,  plastic  straws,  6  cm  long  with  an 
inner  (outer)  diameter  of  0.25  cm  (0.32  cm),  conditioned  the  air  flow  for  even  distribution  over 
the  copper  screen.  Air  flow  rates  of  8  L/min  and  20  L/min,  as  determined  by  a  mass-flow 
controller  (Sierra  Control  Flo-Box  Model  #  905C-PS-BM-11),  were  used  to  generate  foams  with 
expansion  ratios  of  100  and  200,  respectively. 

A  diagram  of  the  drainage  experiment  is  shown  in  Figure  6.  Foam  was  generated 
beneath  the  foam-generation  screen,  moves  through  the  foam  generator,  and  exits  into  a  20- 
liter,  transparent,  plastic  cylinder  with  a  diameter  of  0.213  m  and  a  height  of  0.561  m.  The  base 
of  the  cylindrical  column  was  formed  with  a  sheet  of  aluminum  foil,  sloped  slightly  downward 
toward  the  center,  where  a  1-cm  diameter  hole  allowed  the  drained  liquid  to  pass  out  of  the 
cylinder.  The  liquid  passing  through  the  hole  entered  a  funnel,  collecting  the  drained  liquid  into 
a  500  mL  graduated  cylinder.  A  video  camera  recorded  images  of  the  liquid  filling  the 
graduated  cylinder. 

For  each  foam  drainage  experiment,  the  desired  air  flow  was  first  set  followed  by  the 
liquid.  A  stopwatch  was  used  to  record  the  times  of  key  events.  The  liquid  and  air  combined  on 
the  screen;  the  liquid  foam  solution  was  manually  distributed  on  the  screen  to  start  foam 
generation.  The  excess  liquid,  mostly  from  the  priming  process  before  foam  is  formed,  was 
collected  in  the  passage  beneath  the  screen.  Foam  exited  the  generator  channel  and  traveled 
down  the  side  of  the  cylinder.  The  beginning  of  the  fill  time  was  recorded  when  the  foam 
moved  past  the  specified  height  in  the  collection  cylinder.  Once  the  foam  reached  the  bottom 
of  the  cylinder  and  touched  the  aluminum  foil,  a  video  camera  began  recording  (Pinnacle  Studio 
12)  the  drained  liquid  collecting  in  a  graduated  cylinder.  Once  the  foam  filled  the  20-liter 
cylinder,  the  fill  time  was  recorded,  the  liquid  and  air  flows  were  turned  off,  and  the  foam 
generator  was  set  aside.  A  plastic  lid  was  placed  on  top  of  the  cylinder  to  minimize  foam 
evaporation.  The  excess  liquid  in  the  channel  below  the  foam-generation  screen  was  collected 
and  measured  in  a  graduated  cylinder. 

The  camera  recorded  video  for  about  six  minutes  after  the  foam  reached  the  base  of  the 
cylinder.  The  video  recording  was  used  to  determine  the  amount  of  liquid  drained  versus  time. 
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After  six  minutes,  the  drainage  rate  slowed  significantly,  allowing  the  accumulating  liquid 
volume  in  the  graduated  cylinder  to  be  measured  by  visual  inspection.  The  volume  of  drained 
liquid  was  recorded  at  one-minute  intervals  until  the  change  in  volume  was  less  than  2  mL/min. 
Subsequently,  the  drained  liquid  volume  was  measured  in  intervals  of  several  minutes  until  the 
total  drainage  time  reached  about  an  hour.  The  foam  was  left  in  the  cylinder  overnight  to 
completely  drain,  and  this  liquid  volume  was  used  to  calculate  the  expansion  ratio  of  the  foam. 
Additionally,  a  value  of  the  expansion  ratio  was  also  derived  from  the  net  amount  of  liquid 
sprayed  into  the  foam  during  the  filling  process. 

Volume  of  liquid  in  foam  =  Liquid  flow  rate  x  Fill  time  -  Excess  liquid  (62) 

„  Total  foam  volume 

Ex  = -  (63) 

Volume  of  liquid  in  foam 
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Figure  6:  Experimental  apparatus. 


Foam  bubble  sizes  were  determined  using  photographs  of  the  bubbles  and  a  ruler  taken 
at  the  top  of  the  foam  column  immediately  after  filling.  In  Figure  7,  an  example  photograph 
shows  the  centimeter-size  bubbles  characteristic  of  HiEx  foams. 


29 


Figure  7:  Photograph  of  foam  bubbles  for  Ex  =  290 .  The  ruler  markings  are  in  centimeters. 


The  filling  process  in  the  experiment  contrasts  with  the  idealized  filling  process  of  the 
model,  wherein  the  foam  front  has  a  uniform  height  as  it  travels  upward.  In  the  experiment,  the 
foam  flows  into  the  cylinder  and  a  portion  of  it  slides  along  the  cylinder  wall  as  it  travels 
downward.  After  it  reaches  the  cylinder  base,  the  foam  takes  a  few  seconds  to  reach  the  hole 
where  liquid  can  drain  into  the  graduated  cylinder.  The  liquid  draining  from  the  foam  must 
travel  to  the  hole  in  the  aluminum  foil  and  accumulate  to  form  a  drop,  and  each  drop  must 
travel  down  the  funnel  and  into  the  graduated  cylinder. 

The  experiment  and  model  also  differ  at  the  bottom  boundary.  In  the  model,  liquid 
accumulates  beneath  the  foam.  In  the  experiment,  liquid  drained  from  the  foam  flows  into  a 
graduated  cylinder  beneath  the  foam  column.  The  differences  between  the  model  and 
experiment  are  discussed  below. 

The  liquid  layer  under  the  foam  in  the  model  allows  liquid  to  be  transported  back  into 
the  foam  from  the  water  layer  beneath.  Capillary  suction  is  a  physical  phenomenon  that  occurs 
when  the  pressure  of  the  liquid  in  the  foam  is  less  than  the  pressure  of  the  liquid  with  which  it 
comes  into  contact.  The  suction  of  liquid  can  occur  in  the  model  at  the  beginning  of  the  filling 
process  to  some  extent,  depending  on  the  parameters  of  the  foam.  In  the  experiment, 
however,  there  is  no  liquid  layer  beneath  the  foam.  Under  the  experimental  conditions,  the 
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model  predicts  that  there  only  would  be  a  small  amount  of  liquid  drawn  into  the  foam  over  a 
short  period  of  time  were  the  foam  sitting  on  liquid.  However,  for  much  smaller  bubble  sizes 
and/or  much  greater  expansion  ratios  than  were  experimentally  studied,  the  model  predicts 
that  capillary  suction  is  no  longer  negligible.  Boundary  conditions  have  been  included  in  other 
modeling  works  [8,9]  to  eliminate  capillary  suction  and  introduce  an  induction  period  during 
which  liquid  does  not  drain  from  the  foam.  Because  the  amount  and  duration  of  capillary 
suction  predicted  by  the  model  for  the  experiments  of  this  work  are  very  small,  we  do  not 
correct  the  boundary  conditions  to  eliminate  the  suction. 

The  liquid  drained  during  filling  is  enough  to  form  a  very  small  liquid  layer  beneath  the 
foam  as  it  drains;  thus  the  experimental  bottom  boundary  condition  during  filling  should  be 
similar  to  that  of  the  model.  However,  the  experimental  foam  may  not  have  a  thin  layer  of 
liquid  present  during  the  free  drainage  period.  As  a  result,  the  model  could  slightly  overpredict 
the  effect  of  capillarity  during  the  free  drainage  process.  Additionally,  the  model  considers  the 
reduction  of  the  length  of  the  foam  layer  due  to  the  increasing  height  of  the  water  layer.  The 
corresponding  reduction  in  height  is  negligible  for  HiEx  foams  because  the  height  of  the  water 
layer  is  very  small  compared  to  the  height  of  the  foam  layer,  i.e.  hw  «  hf .  This  is  an 

approximation  made  in  the  analytical  solutions  of  the  model. 

Finally,  the  foams  generated  in  experiments  have  a  distribution  of  bubble  radii,  whereas 
the  model  assumes  all  bubbles  are  equal  in  size.  A  direct  comparison  between  polydisperse  and 
monodisperse  foams  is  not  possible,  and  it  is  not  clear  what  averaging  procedure  will  provide  a 
mean  bubble  radius  that  will  correlate  with  the  experimental  drainage  data  for  the  polydisperse 
foam.  Magrabi  et  al.  [8]  propose  that  the  appropriate  average  bubble  size  is  calculated  from  the 
bubble-size  dependence  of  the  liquid  flow  rate  through  a  channel.  In  their  work,  the  average 
liquid  flow  rate  through  a  channel  in  a  polydisperse  foam  is  calculated  by  the  fourth  moment  of 
the  bubble  radius  distribution  because  the  flow  rate  is  proportional  to  the  fourth  power  of  the 
bubble  radius.  In  our  experiments,  we  did  not  obtain  a  reliable  bubble-size  distribution  function 
that  would  enable  us  to  calculate  an  average  bubble  radius.  Instead,  we  adjusted  the  bubble 
radius  to  fit  the  experimental  data  for  each  foam,  and  the  bubble  radius  obtained  from  the  fit 
was  compared  to  the  range  of  bubble  radii  obtained  from  photographs  such  as  that  shown  in 
Figure  7.  For  all  experiments,  the  bubble  radius  obtained  from  the  fit  was  within  the  range  of 
bubble  radii  observed  in  the  photographs. 


4.0  RESULTS  AND  DISCUSSION 

An  example  case  of  foam  drainage  during  filling  is  shown  in  Figure  8.  At  the  start  of 
filling,  there  is  a  relatively  large  pressure  gradient  caused  by  the  small  foam  height  and 
boundary  conditions  on  a  that  competes  with  gravity.  The  drainage  behavior  during  this  initial 
period  strongly  depends  on  the  choice  of  parameters.  As  the  foam  height  increases,  the  force  of 
gravity  becomes  dominant  and  drives  liquid  out  of  the  foam.  A  constant  drainage  rate  from  the 
foam  quickly  develops,  and  the  height  of  the  water  layer  increases  linearly  with  time.  During 
this  period,  the  average  liquid  fraction  of  the  foam  decreases  to  a  steady-state  value.  For  the 
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example  case  shown,  the  average  liquid  fraction  decreases  from  an  input  value  of  0.005 
(Ex  =  200)  to  a  steady-state  value  of  about  0.00425  (Ex  =  235)  .  The  foam  loses  about  15%  of 
its  original  liquid  during  the  filling  process. 

Under  other  conditions,  such  as  slower  injection  velocities,  very  high  expansion  ratios, 
and/or  small  bubble  sizes,  hw  decreases  from  zero  at  the  beginning  of  the  filling  process.  During 
this  time,  the  pressure  gradient  dominates  gravity  and  causes  liquid  suction  into  the  foam  from 
the  water  layer  beneath.  This  behavior  is  also  reflected  in  the  behavior  of  a ,  which  rises  above 
its  initial  value  at  the  start  of  filling. 
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Figure  8:  Model  prediction  of  liquid  drained  during  foam  filling. 
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Figure  9:  Model  prediction  of  liquid  drained  during  and  after  foam  filling. 


At  the  end  of  filling,  the  free  drainage  period  begins.  Figure  9  extends  the  computation 
in  Figure  8  to  include  the  free  drainage  behavior  after  filling.  The  foam  continues  to  drain,  but  is 
no  longer  receiving  liquid  from  input  foam.  When  filling  is  ceased,  the  drainage  rate  begins  to 
decrease  slightly  and  the  average  liquid  fraction  decreases  abruptly.  In  the  foam,  the  channels 
contract  in  response  to  the  reduced  liquid  content,  and  the  drainage  rate  will  approaches  zero 
as  the  driving  forces  for  liquid  flow  come  into  equilibrium.  During  this  process,  the  height  of  the 
water  layer  approaches  its  equilibrium  value  at  an  exponential  rate,  and  the  average  liquid 
fraction  behaves  in  a  similar  fashion. 

The  analytical  solutions  of  this  work  are  valid  only  for  a  subset  of  the  conditions  for 
which  their  original  differential  equations  can  be  solved  by  numerical  methods.  For  the  filling 
case,  the  approximations  made  in  the  analytical  solution  are  valid  when  the  average  liquid 
fraction  decreases  from  the  input  liquid  fraction  by  much  less  than  50%.  It  is  difficult  to  predict 
a  priori  when  this  condition  will  be  satisfied.  To  show  the  important  parameters  to  consider 
when  deciding  whether  to  use  the  analytical  solution,  we  consider  Eq.  (30)  when  hw  «  h,  and 

the  linear  term  in  Eq.  (43)  to  yield  the  approximation 
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From  the  previous  expression,  it  can  be  seen  that  the  quantity  affecting  the  reduction  in  a 
from  the  input  liquid  fraction  is  vinj  /  Aaf .  By  comparison  of  the  drainage  rate  as  predicted  by 

the  numerical  and  analytical  solutions,  we  see  that  the  difference  for  Ex  up  to  200  is  around 
2%  when  we  use  the  condition  vinj  /  Aaf  ~  5  .  Note  that  for  Ex  up  to  600,  a  slightly  greater  value 

of  8  is  needed  to  achieve  the  same  accuracy.  Substitution  of  this  value  in  Eq.  (64)  indicates  that 
the  average  liquid  fraction  should  not  decrease  from  the  input  value  by  much  more  than  about 
15%  for  the  analytical  solution  to  maintain  this  accuracy.  For  our  experiments,  the  value  of 
vmj  /  Aaf  ranges  from  about  0.5  to  2,  and  the  analytical  solution  does  not  provide  the  drainage 

rate  with  sufficient  accuracy. 

It  is  important  to  note  that  the  analytical  expression  is  valid  for  field  cases,  where  the  fill 
rate  is  much  faster.  For  a  foam  delivered  to  a  room  with  a  standard  ceiling  height  of  H  «  2.5  m 
that  must  be  filled  within  120  s ,  one  obtains  v,„,  /  Aa,  ~  10  when  using  conditions  similar  to 
our  experimental  values,  A  «  0.45  ml s  and  af=  1  / 200  .  In  Figure  10,  we  show  how  the 
analytical  solution  compares  with  the  numerical  solution  for  bench-scale  and  large-scale  values 
°f  Viv  /  Aaf . 
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Figure  10:  Comparison  of  numerical  and  analytical  solutions  for  the  height  of  the  liquid  layer 
as  a  function  of  time  during  filling  for  large-scale  and  bench-scale  conditions.  The  parameters 
for  the  large  scale  are  v.  .  =  0.05  mis  and  af  =  0.002 .  The  large-scale  injection  velocity  is 

calculated  based  on  a  6  m  column  height  and  a  120  s  fill  time.  The  parameters  for  the  bench 
scale  are  vinJ  =  0.0022  ml s  and  af  =  0.0083  .  In  both  calculations,  the  parameters  £  =  0.0066 

and  Db  =  0.006  m  are  used. 


The  analytical  solution  for  free  drainage  is  in  excellent  agreement  with  the  numerical 
solution  to  the  original  differential  equation.  In  Figure  11,  we  show  the  drainage  behavior  for 
different  bubble  sizes  in  the  case  of  initial  Ex  =  200 .  the  solutions  are  nearly  identical.  Initial 

expansion  ratio  values  were  also  varied  independently  of  bubble  size  from  Ex  =  20  to  600  and 
the  same  excellent  agreement  was  obtained. 
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Figure  11:  Comparison  of  numerical  and  analytical  solutions  for  the  height  of  the  drained 
liquid  for  free  drainage  for  different  bubble  sizes. 


Shown  in  Figure  12  is  a  comparison  of  solutions  for  filling  and  free  drainage  combined. 
The  analytical  solutions  using  the  large-scale  parameters  are  identical  to  the  numerical 
solutions.  However,  the  analytical  solutions  calculated  with  the  bench-scale  parameters  exhibit 
significant  differences  during  both  filling  and  free  drainage  periods.  The  overprediction  for 
drainage  during  filling  is  responsible  for  the  underprediction  in  the  free  drainage  period 
compared  to  the  numerical  solution.  The  initial  conditions  for  hw  and  a  in  the  analytical  free- 
drainage  solution  are  equal  to  the  corresponding  values  of  the  filling  analytical  solution  at  the 
end  of  the  filling  process.  Hence,  it  is  necessary  to  accurately  calculate  the  filling  solution  to  also 
accurately  predict  the  magnitude  of  the  liquid  drained  during  the  free-drainage  period. 

The  model  indicates  that  the  drainage  rate  from  the  filling  solution  approaches  the 
initial  drainage  rate  of  the  free  drainage  solution  at  large  values  of  v. .  /  Aaf .  The  difference 

between  the  initial  drainage  rates  at  vinj  /  Aaf  =20  is  about  10%,  and  at  vinj  /  Aaf  =40,  the 

difference  is  about  5%.  This  suggests  that  filling  can  be  ignored  at  very  high  injection  velocities 
because  the  analytical  free-drainage  solution  provides  a  reasonable  approximation  to  the  entire 
drainage  process. 
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Figure  12:  Comparison  of  numerical  and  analytical  solution  both  during  and  after  the  filling 
process.  The  parameters  for  the  large  scale  are  vinJ  =  0.05  mis,  tjm  =  120  s ,  H  =  6  m ,  and 

af  =  0.002 .  The  parameters  for  the  bench  scale  are  v  =  0.0022  mis,  tfUl  =  259  s , 

H  =  0.561m ,  and  af  =  0.0083 .  In  both  calculations,  the  parameters  E,  =  0.0066  and 
Db  =  0.006  m  are  used. 


The  analytical  solutions  of  this  work,  containing  many  approximations,  can  be  used  to 
calculate  drainage  behavior  in  good  agreement  with  numerical  solutions  of  the  partial 
differential  equations  from  which  the  model  is  derived.  To  show  this,  we  compare  our  results 
with  those  of  another  model  that  has  been  used  to  describe  free  drainage,  that  of  Magrabi  et 
al.  [8],  which  uses  similar  starting  equations. 

Magrabi  et  al.  [8]  studied  free  drainage  from  foam  with  an  initial  expansion  ratio  of  20 
(a  =  0.05) .  In  Figure  13,  we  compare  their  modeling  results  with  modeling  results  using  our 
simplified  analytical  solution  for  free  drainage.  The  parameters  are  listed  in  Appendix  B. 

For  comparison,  we  use  the  data  from  Ref.  [8]  for  the  mass  of  water  drained  as  a 
function  of  time  to  calculate  the  corresponding  value  of  hw,  the  volume  of  liquid  drained  per 
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unit  area  of  the  foam  column  used  in  the  experimental  work.  Using  the  liquid  density  and  foam- 
column  geometry  reported  in  Ref.  [8],  we  transform  the  drained-mass  data  to  water  height  via 


K(?)  = 


Ml 

s 


Md{t ) 

pnRj 


(65) 


where  Vd(t )  is  the  volume  of  liquid  drained  from  the  foam,  S  is  the  cross-sectional  area  of  the 
foam  column,  M  d(t )  is  the  mass  of  liquid  drained,  p  is  the  intrinsic  liquid  density,  and  R  is  the 
radius  of  the  cylindrical  foam  column. 

We  shifted  by  the  respective  induction  times  as  reported  in  Ref.  [8]  of  Magrabi  et  al.  In 
the  work,  there  was  a  no-flow  boundary  condition  at  the  bottom  of  the  foam  column  until  that 
boundary  reached  a  liquid  fraction  of  0.26.  Our  model,  however,  has  a  constant  bottom 
boundary  condition  of  aw  =  0.26  .  The  induction  time  predicted  by  our  model  is  much  smaller, 
in  comparison. 


Figure  13:  Height  of  drained  liquid  for  free  drainage  (experimental  and  modeling  results  from 
Ref.  [8]  and  modeling  results  of  this  work).  The  induction  time  of  Ref.  [8]  is  tind  =  240  s  .for 
Dh  =  266  pm  and  tmd  «  20  s  for  Dh  =  540  pm .  The  induction  time  for  our  modeling  results  is 
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As  shown  in  Figure  13,  the  experimental  data  of  Ref.  [8]  disagree  with  the  predictions  of 
their  model.  The  authors  argue  that  coarsening  is  the  effect  responsible  for  the  disagreement, 
which  is  supported  by  measurements  of  average  bubble  radius  that  increase  significantly  during 
the  first  960  s  after  foam  formation.  The  effect  of  coarsening  was  not  considered  in  their 
model,  but  it  was  shown  that  by  increasing  the  bubble  size  from  the  initial  value  to  the  value 
measured  after  960  s,  the  model  predicts  a  greater  drainage  rate  that  increases  the  agreement 
with  the  experimental  data.  However,  bubble  sizes  are  difficult  to  measure  during  drainage 
because  the  distribution  evolves  over  time.  Further,  the  physics  responsible  for  the  evolution  of 
bubble  size  in  foam  is  not  well  understood. 

Although  our  model  calculations  of  their  conditions  show  disagreement  with  their 
experimental  data,  our  modeling  calculations  are  in  reasonable  agreement.  The  experimental 
foam  in  Ref.  [8]  contained  a  low  expansion  ratio  of  20  and  a  bubble  diameter  less  than  a 
millimeter;  our  model  was  designed  specifically  to  treat  high-expansion  foams  (Ex  >  100)  with 
bubble  diameters  ranging  from  several  millimeters  to  a  few  centimeters.  Coarsening  will  have 
much  less  of  an  effect  on  the  drainage  behavior  for  our  foams  with  bubble  diameters  an  order 
of  magnitude  larger.  Additionally,  our  model  predicts  a  smaller  drainage  rate  than  their  model. 
The  largest  contribution  to  the  difference  in  the  model  predictions  is  that  we  consider  the 
length-averaged  behavior  of  the  liquid  volume  fraction  and  liquid  velocity  by  integrating  the 
conservation  equations.  The  resulting  simplification  of  the  model  was  expected  to  introduce 
some  error.  However,  the  clear  advantage  of  considering  the  average  behavior  is  that  it  allows 
the  calculation  of  analytical  solutions  to  describe  the  drainage  behavior  as  a  function  of  time. 

The  liquid  fraction  predictions  of  the  two  models  are  compared  in  Figure  14.  The  figure 
shows  how  the  height-averaged  liquid  fraction  predicted  in  this  work  compares  to  the 
predictions  of  a  height-resolved  model  at  various  heights  within  a  column.  From  visual 
inspection,  our  analytical  solution  appears  to  capture  the  length-averaged  behavior  reasonably 
well. 
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Figure  14:  Comparison  of  length-averaged  liquid  fraction  a  of  this  work  with  model 
predictions  of  liquid  fractions  at  four  positions  along  the  height  of  the  column  from  Ref.  [8]. 
Note  that  z/H=0  is  the  bottom  of  the  column  and  z/H=l  is  at  the  top. 


Next,  we  provide  a  comparison  of  our  model  with  experimental  drainage  measurements 
during  both  filling  and  free  drainage.  We  reiterate  that  the  analytical  solution  for  filling  is  not 
valid  for  the  relatively  slow  injection  velocities  in  the  bench-scale  experiments  of  this  work, 
although  it  is  valid  for  large-scale  field  conditions.  Therefore,  during  the  filling  portion  of  the 
experiments,  we  numerically  solve  Eq.  (31)  to  obtain  hw  as  a  function  of  time  from  t  -  0  until 
the  fill  time.  On  the  other  hand,  the  analytical  solution  for  free  drainage  is  valid  for  the  bench- 
scale  experiments,  and  we  use  the  analytical  solutions  for  1  greater  than  the  fill  time. 
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Figure  15:  Comparison  of  model  and  experiment  during  the  foam-filling  period.  The  induction 
time  is  denoted  tind  . 


The  modeling  and  experimental  data  of  this  study  for  filling  are  shown  in  Figure  15.  The 
experimental  data  suggest  that  the  height  of  the  water  layer  increases  linearly  with  time,  as 
predicted  by  the  model.  Note  that  there  is  an  induction  time  of  approximately  20  s  in  the 
experiments  of  this  work,  but  the  model  underpredicts  this  value. 

Note  the  appreciable  drainage  that  occurs  during  filling  in  Figure  15.  For  the  Ex=183 
case,  the  length-averaged  liquid  fraction  decreases  from  about  0.0055  (Ex  =  183)  to 
approximately  0.00425  (Ex  =  235)  during  the  filling  period  of  103  s.  In  the  Ex  =  120  case,  the 
liquid  fraction  decreases  from  0.0083  (.fix  =  120)  to  about  0.0046  (Ex-  217)  during  the  filling 
period  of  259  s.  The  expansion  ratio  increases  by  about  28%  and  81%,  respectively,  during  the 
filling  period.  This  shows  that  the  drainage  occurring  during  the  filling  process  can  significantly 
affect  the  expansion  ratio  of  the  foam. 

The  model  contains  the  adjustable  parameter  £ .  We  assume  that  the  HiEx  foams  of  this 
study  have  immobile  channel  walls  and  we  use  a  permeability  constant  for  these  conditions 
that  has  been  calculated  in  Ref.  [7],  Also,  the  bubble  diameter  for  the  model  calculations  is  fit  to 
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the  experimental  data  of  liquid  drained  versus  time.  We  fit  the  filling  portion  of  the 
experimental  drainage  curve  because  we  use  a  modified  Hagen-Poiseuille  equation,  valid  for 
fully  developed,  steady-state  flow.  There  is  a  wide  distribution  in  the  experimental  bubble  size, 
ranging  in  diameter  from  about  a  millimeter  to  a  couple  of  centimeters.  The  average  bubble 
size  from  fitting  the  experimental  data  was  shown  to  be  within  the  bubble-size  distribution 
from  the  photographs.  The  effect  of  bubble  size  on  drainage  behavior  is  examined  below  and 
the  bubble-size  range  within  the  distribution  observed  in  the  experiments  is  investigated. 
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Figure  16:  Comparison  of  model  with  experimental  data  for  Ex=  183  for  filling  and  free 
drainage  periods. 


The  experimental  data  and  model  predictions  for  the  entire  drainage  period  including 
both  filling  and  free  drainage  are  shown  in  Figure  16  and  Figure  17  for  the  cases  of  Ex  - 183 
and  Ex  =  120 ,  respectively.  After  the  filling  period,  the  drainage  rate  in  the  experiment 
decreases.  Good  agreement  is  shown  between  the  model  and  experiment  during  the  free 
drainage  period,  which  suggests  that  the  data  fits  well  to  the  exponential  dependence  of  the 
model. 
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Figure  17:  Comparison  of  model  and  experimental  data  for  Ex=120  during  filling  and  free 
drainage  periods. 


Comparison  of  model  and  experiment  results  of  this  work  show  better  agreement  than 
the  model  and  experiment  results  comparison  of  Magrabi  et  al  [8],  In  the  latter,  the  predicted 
value  of  hw  at  equilibrium  is  50%  smaller  than  the  experimental  value.  Our  predicted  value  of 
hw  at  equilibrium  is  very  close  to  the  experimental  value.  To  understand  the  difference,  we 
compare  the  ratio  of  the  equilibrium  average  liquid  fraction  to  the  initial  liquid  fraction  for  the 
two  experiments.  We  use  Eq.  (60)  to  calculate 
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Note  that  this  is  also  equal  to  the  ratio  of  the  liquid  volume  of  the  foam  at  equilibrium  to  the 
initial  liquid  volume  in  the  foam.  Using  the  parameters  of  Magrabi  et  al.  [8]  for  a  bubble  radius 
of  133  pm  and  JaJ  «  y[a~, ,  this  ratio  is  approximately  0.5;  only  50%  of  the  water  is  drained 

from  the  foam.  Note  that  the  equilibrium  experimental  value  in  Ref.  [8]  is  approximately  equal 
to  a(0)H  =  0.05  x  0.2  in  —  0.01  m ,  the  total  liquid  volume  in  the  foam  per  unit  area  of  the  foam 


43 


column,  and  therefore  the  experiment  indicates  that  nearly  100%  of  the  liquid  in  the  foam  is 
drained  at  equilibrium.  The  prediction  of  the  percent  of  liquid  drained  at  equilibrium  using  Eq, 
(66)  is  consistent  with  the  model  data  shown  in  Figure  13;  the  value  of  hw  is  about  50%  of  the 
experimental  value.  On  the  other  hand,  this  ratio  is  approximately  0.02  for  the  parameters  of 
this  study;  only  2%  of  the  initial  liquid  remains  in  the  foam.  The  main  contribution  to  the 
difference  in  the  a  values  is  the  very  large  difference  in  bubble  size  between  the  two  studies. 


Table  2:  Model  parameters  for  comparison  with  experimental  data. 


r 

Surface  tension 

7  N 

2.25  xlCT2  — 
m 

A 

Liquid  viscosity 

8.94 xKT4  kg 
m  s 

4 

Permeability  constant 

0.0051 

p 

Liquid  density 

997 

m 

Liquid  fraction  at  foam/liquid  interface 
(mass  balance) 

1 

«vv 

Liquid  fraction  at  foam/liquid  interface 
(momentum  balance) 

0.26 

Table  3:  Initial  conditions  for  filling  solution 


K o 

Initial  height  of  liquid  layer 

0  m 

hfo 

Initial  height  of  foam  layer 

0.01  m 

*0 

Starting  time 

0  s 

An  advantage  of  the  analytical  solutions  of  this  work  is  that  they  can  be  used  to  easily  fit 
the  permeability  coefficient  to  experimental  data  for  a  given  foam.  In  this  work,  the 
permeability  coefficient  appears  to  be  constant  between  two  experiments  where  the  expansion 
ratio  is  varied  from  183  to  120  and  the  fill  time  is  varied  from  103  s  to  259  s  in  a  foam  column  of 
height  0.561  m.  Further  work  is  necessary  to  determine  the  range  in  foam  parameters  for  which 
the  permeability  coefficient  remains  approximately  constant.  Note  that  previous  models  [6,7] 
which  assume  that  the  viscous  force  is  predominantly  located  in  channels  with  immobile  walls 
also  have  an  effective  permeability  coefficient  that  is  constant,  regardless  of  foam  parameters. 
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Figure  18:  Effect  of  fill  rate  on  drainage. 


EFFECT  OF  MODEL  PARAMETERS  ON  FOAM  DRAINAGE  BEHAVIOR 
A.  Fill  Rate 

The  effect  of  fill  rate  on  drainage  is  shown  in  Figure  18.  A  case  with  a  very  small  fill  time 
of  5  s  is  compared  with  cases  with  slower  injection  velocities  corresponding  to  fill  times  of  120  s 
and  600  s  (10  minutes).  For  the  height  of  foam  column  and  parameters  considered,  the 
drainage-time  difference  between  filling  in  5  s  and  120  s  is  small.  However,  spaces  for  which 
HiEx  foams  are  typically  used  for  fire  suppression  have  higher  ceilings  and  thus  the  foam 
column  heights  are  potentially  taller  by  an  order  of  magnitude.  Hence,  the  fill  times  could  be 
much  larger.  The  curve  corresponding  to  a  600s  fill  time  shows  a  clear  difference  from  the  other 
curves;  this  difference  cannot  be  predicted  using  a  free-drainage  model. 
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Figure  19:  Effect  of  bubble  size  on  drainage. 


B.  Bubble  Size 

The  effect  of  the  bubble  size  on  drainage  is  shown  in  Figure  19.  The  model  predicts  that 
foams  with  larger  bubbles  drain  at  a  faster  rate  and  drain  more  liquid  before  reaching 
equilibrium.  This  effect  has  been  observed  in  previous  studies  at  lower  expansion  ratios  near  20 
and  bubble  sizes  less  than  a  millimeter  [9],  At  a  given  liquid  fraction,  a  foam  with  larger  bubbles 
will  have  channels  with  greater  cross-sectional  area.  Because  the  average  velocity  in  the 
channel  is  proportional  to  the  channel  area,  the  drainage  rate  from  the  foam  will  be  greater. 

In  Figure  19,  drainage  is  predicted  for  a  bubble  diameter  of  0.001  m  (1  mm).  In  studies 
of  low-expansion  foams,  the  bubble  diameters  are  typically  smaller.  The  model  indicates  that 
Hi  Ex  foams  do  not  drain  appreciably  for  these  smaller  bubble  sizes.  Flowever,  the  foams  used  in 
practice  for  HiEx  can  have  bubble  diameters  up  to  a  few  centimeters.  As  Figure  19  indicates. 

Hi  Ex  foams  drain  appreciably  during  and  after  the  filling  process  for  the  bubble  sizes  used  in 
practice. 
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Figure  20:  Effect  of  input  expansion  ratio  on  drainage. 


C.  Expansion  Ratio 

The  effect  of  varying  the  expansion  ratio  of  the  input  foam  on  drainage  behavior  is 
displayed  in  Figure  20.  For  a  given  bubble  size,  the  area  of  the  channels  will  decrease  with 
increasing  expansion  ratio.  The  drainage  rate  is  therefore  reduced  in  foams  with  higher 
expansion  ratios. 
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Figure  21:  Effect  of  column  height  on  drainage. 


D.  Column  Height 

The  effect  of  column  height  on  drainage  is  shown  in  Figure  21.  As  the  height  of  the  foam 
column  increases,  the  time  required  to  reach  equilibrium  also  increases.  Further,  the  liquid 
fraction  at  equilibrium  also  decreases  with  increasing  column  height.  This  behavior  is  caused  by 
the  pressure  gradient  of  the  model,  which  decreases  with  increasing  column  height.  The  effect 
of  capillarity  on  the  column  will  therefore  be  smaller  for  taller  columns,  allowing  a  greater 
percentage  of  liquid  to  drain  before  equilibrium  is  reached. 
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Figure  22:  Effect  of  permeability  coefficient  on  drainage. 


E.  Permeability  Coefficient 

The  effect  of  varying  the  permeability  coefficient  on  drainage  is  plotted  in  Figure  22.  The 
drainage  behavior  for  the  permeability  constant  E,  =  0.0066  is  shown  along  with  greater  values, 
0.013  and  0.030.  These  values  correspond  approximately  to  conditions  that  have  been  reported 
in  the  literature.  Koehler  et  al.  [7]  reported  a  calculated  value  of  a  permeability  constant  for 
channels  with  immobile  walls,  E,  «  0.0063 ,  which  is  the  similar  in  value  to  the  permeability 
constant  used  in  our  work.  Magrabi  et  al.  [9]  used  a  permeability  coefficient  of  about  0.025  and 
0.01  to  model  foams  at  an  expansion  ratio  of  20  made  from  the  surfactant  solutions  aqueous 
film-forming  foam  (AFFF)  and  film-forming  fluoroprotein  foam  (FFFP),  respectively. 


5.0  CONCLUSIONS 

We  have  developed  a  model  to  describe  the  loss  of  liquid  from  foam  that  includes  an 
essential  process  of  industrial  foam  applications:  the  process  of  filling  a  space  with  foam.  The 
model  solves  for  the  time  evolution  of  liquid  drained  from  the  foam  and  the  average  liquid 
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content  within  the  foam.  Introducing  the  filling  process  adds  a  new  time  scale,  the  fill  time,  to 
the  drainage  problem  that  has  a  significant  effect  on  the  drainage  behavior. 


An  important  aspect  of  the  work  is  that  approximations  in  the  model  equations  that 
describe  the  drainage  of  liquid  from  foam  are  made  to  provide  simple  analytical  solutions.  The 
analytical  solutions  indicate  that  the  drainage  rate  remains  constant  during  the  filling  process 
and,  after  filling,  decays  exponentially. 

The  investigation  is  aimed  at  describing  drainage  from  foams  with  very  low  liquid 
content,  such  as  those  used  in  high-expansion  (HiEx)  foam  fire-suppression  applications.  Model 
results  are  in  good  agreement  with  the  experimental  data  from  bench-scale  experiments  on 
HiEx  foams.  Although  numerical  solutions  were  needed  for  the  model  equations  to  obtain 
sufficient  accuracy  under  the  bench-scale  conditions,  the  simple  analytical  solutions  of  the 
model  are  shown  to  be  valid  for  field  conditions. 

Modeling  predictions  were  compared  with  experimental  data  using  a  permeability 
coefficient,  £  =  5.1xlCT3,  that  was  calculated  from  the  permeability  derived  in  Ref.  [7]  for  a  no¬ 
slip  condition  at  the  channel  walls.  The  model  was  fit  to  experimental  data  by  adjusting  the 
bubble  size.  The  bubble  size  that  provided  the  best  fit  to  experiment  results  was  within  the 
bubble-size  distribution  obtained  by  photographs  of  the  experimental  foams.  Fill  rate,  bubble 
size,  input-foam  expansion  ratio,  column  height,  and  permeability  coefficient  each  have  a 
substantial  impact  on  the  foam  drainage  behavior  predicted  by  the  model. 
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APPENDIX  A:  Calculation  of  the  superficial  liquid  velocity  for  a  cylindrical  channel 

In  order  to  treat  liquid  flow  in  foam,  we  first  assume  the  foam  is  composed  of  a  network 
of  cylindrical  channels.  We  calculate  the  superficial  liquid  velocity  for  a  cylindrical  channel  to 
show  its  dependence  on  the  channel  area.  In  the  calculation,  we  consider  steady-state  flow. 
Also,  we  assume  that  the  walls  are  immobile,  i.e.  there  is  a  no-slip  condition  at  the  surface  of 
the  channels.  This  gives  rise  to  a  parabolic  velocity  profile.  The  velocity  of  the  liquid  within  the 
channel  at  a  distance  r  from  the  axis  can  be  described  by  the  equation 


(67) 


where  rp  is  the  radius  of  the  channel. 


The  stress  at  a  distance  r  is  given  by 
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(68) 


We  are  interested  in  modeling  the  stress  at  the  surface  of  the  channel  as  an  interaction 
between  phases  in  the  multiphase  momentum  equation.  Hence,  we  want  the  value  of  the 
stress  at  r  =  rp 


Trz (rp )  =  .  (69) 

rr 

We  must  calculate  the  force  from  this  stress  by  integrating  over  the  surface  area  of  the  channel 

F„  =  J  dA  =  \}nrpL\  =  4 n/uvmaxL  (70) 

rp 

where  2 nr  L  is  the  surface  area  of  the  channel.  Eq.  (70)  gives  the  force  in  a  single  channel.  In 

order  to  find  the  total  force  per  unit  volume  of  the  simplified  foam,  we  need  to  multiply  by  the 
number  of  channels  per  unit  volume. 

For  dodecahedral  geometry,  there  are  np  =  10  channels  per  bubble.  For  this  case,  the  frictional 
force  per  unit  volume,  Ff ,  is 


Ff=FrznpNh=lOFr:Nh 


(71) 
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where  Nb  is  the  number  of  bubbles  per  unit  volume.  We  can  calculate  the  total  number  of 
bubbles,  Nb ,  by 


Vf 

(72) 

"b 

where  Vf  is  the  total  volume  of  the  foam  and  Vb  is  the  volume  per  bubble.  Dividing  the 
previous  expression  by  V r  yields 
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The  force  per  unit  volume  due  to  friction  at  the  surface  of  the  channels  is  therefore 
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We  want  to  model  this  frictional  force  as  an  interaction  between  phases  in  the  liquid-phase 
momentum  equation.  The  momentum  equation,  ignoring  the  inertial  and  accumulation  terms, 
for  the  liquid  phase  (denoted  by  / )  within  a  two-phase  mixture  including  gas  (denoted  by  a  )  is 
given  by 

0  =  -VP0)  -[v-r]</>  +  p(,)g  +  pU)F,a  (va  -v,)  (75) 


where  pU)  =  ap  is  the  effective  density  of  the  liquid  within  the  mixture  and  P(,)  =  aP:  is  the 
effective  liquid  pressure.  The  last  term  in  the  previous  expression  is  the  interaction  between 
the  gas  and  liquid  phases,  and  this  is  the  term  for  which  we  want  to  model  the  friction  between 
liquid  and  gas  in  a  foam  using/y  =  p(,)F,a  (va  -  v,).  The  quantity  Fla  is  termed  a  time  constant 
for  the  interaction  between  the  phases.  The  interaction  term  has  the  form 


-apF\ 


1  (F^pvmax 
(l5  +  7V5)l2 


(76) 


when  we  consider  the  gas  velocity  va  =  0  .  We  also  express  the  liquid  velocity  in  terms  of  the 
superficial  velocity,  v^.To  relate  vmax  to  vs,  let  us  consider  the  pore  velocity,  vp ,  which  is  the 
average  liquid  velocity  within  a  channel.  The  average  velocity,  vp ,  is  related  to  the  maximum 
velocity,  vmax ,  by  a  constant;  let  us  make  the  relation 

v  —cv  .  (77) 

max  p  '  ' 
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We  are  considering  steady-state  flow  in  cylindrical  channels,  and  therefore  c  =  2  for  the 
resulting  parabolic  velocity  profile.  We  also  must  relate  the  channel  velocity  to  the  superficial 
velocity.  We  achieve  this  by  means  of  a  mass  balance  in  a  foam  column  at  the  interface 
between  a  foam  and  liquid  layers.  Let  us  denote  N  as  the  total  number  of  channels  that  are 

connected  and  drain  directly  to  the  foam  layer,  each  having  a  flow  rate  of  apvp  .  The  liquid 
drained  from  the  foam  causes  the  liquid  layer  to  rise,  and  the  flow  rate  of  the  liquid  is  vsA, 
where  A  is  the  cross-sectional  area  of  the  foam  column.  Balancing  these  flow  rates,  we  have 
the  relation 


NpapVp  =  VsA 


VP  = 


vNpaPJ 


(78) 


A  common  equation  in  porous  media  is  vs  =  evp,  where  e  is  the  porosity  of  the  medium.  If  we 

consider  foam  as  a  porous  media  with  liquid  filling  the  pores,  then  e  =  a .  This  allows  the 
relation  a  =  N  a  /  A  .  The  time  constant  is  therefore 
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Our  momentum  equation  for  the  liquid  phase  is 
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Solving  for  vs  yields 


v.  = 


(l5  +  7V5)l2 
160  npc 


A 


NpQp  1  \pag-VP(l) )  . 


If  we  set  N  a  /  A  =  a ,  then 
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The  superficial  velocity  in  terms  of  Darcy's  law  is  given  by 
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(81) 
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Assuming  that  V(aP,)l a  ~  VP, ,  the  previous  expression  can  be  used  to  find  the  permeability  of 
the  foam  with  cylindrical  cross-section  using  Eq.  (82),  yielding 


k  = 


'(15  +  7V5 

160  k/jc 

v  7 


a1!} 


(84) 


This  shows  that  the  permeability  (and  the  superficial  velocity)  is  proportional  to  the  square  of 
a  .  The  permeability  constant  for  this  case  is 


4  = 


(l5  +  7V5) 
160;rc 

V  J 


0.03 


(85) 


The  case  considered  here  is  for  vertically  oriented  cylindrical  channels  with  immobile  walls.  By 
considering  the  equation  for  Hagen-Poiseuille  flow  in  a  cylindrical  channel 


v  =  — - - (cylindrical)  (86) 

8;r  dz 

and  the  equation  for  flow  in  a  channel  with  equilateral-triangular  cross  section  derived  by  Desai 
and  Kumar 


aP  dP 
20^3  dz 


(equilateral  triangular) 


(87) 


we  calculate  a  correction  factor,  /  20>/3  «  0.726 ,  that  can  correct  the  value  in  Eq.  (85)  for 
vertically  oriented  cylindrical  channels  to  yield  £  for  vertically  oriented  channels  with 
equilateral  triangular  cross  section: 

E,  ~  0.03x0.726  w  0.022  (equilateral  triangular)  .  (88) 

The  channels  in  foam  have  a  cross-sectional  area  that  is  shaped  like  a  triangle  wherein  the  three 
edges  are  curved  in  a  concave  manner,  i.e.  a  scalloped-triangular  shape.  We  expect  that  the 
value  of  £  for  vertically  oriented  channels  with  scalloped-triangular  cross  section  and  immobile 
walls  is  slightly  less  than  the  value  provided  in  Eq.  (88). 
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APPENDIX  B:  Parameters  for  comparison  with  Reference  [8], 

The  following  parameters  were  used  to  compare  the  model  in  this  work  with  the  model 
calculations  reported  in  Ref.  [8]  as  shown  in  Figure  13  and  Figure  14. 


Y 

Surface  tension 

,  N 

2.25  xl(T2  — 
m 

P 

Liquid  viscosity 

8.94x10  4  kg 
m  s 

4 

Permeability  constant 

0.0245 

L 

Plateau  border  length 

Note:  L  =  0.408Dfe 

1.09x10 

[Db=  2.66x10  4  m) 

2.20x10^  m 
(Dh  =  5.40x10  4  m) 

P 

Liquid  density 

997 

m 

H 

Column  height 

0.2  m 

«/(=  a  0) 

Liquid  fraction  of  input  foam 

0.05 

«v. 

Liquid  fraction  at  foam/liquid 
interface  (mass  balance) 

1 

Liquid  fraction  at  foam/liquid 
interface  (momentum  balance) 

0.26 

Starting  time 

05 
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